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ABSTRACT 


The  extreme  speed  and  altitude  ranges  of  which  modern,  high  perfor¬ 
mance,  supersonic  airplanes  are  capable  have  become  manifested  in  wide 
variations  in  airframe  transfer  function  parameters  as  flight  conditions 
change*  Airframe  control  systems  and  stability  augmentation  devices 
utilized  in  these  aircraft  require  some  form  of  gain  changing  to  pro¬ 
vide  satisfactory  stability  and  control  characteristic  over  the  complete 
flight  envelope.  In  present  day  operational  aircraft,  system  gains  are 
scheduled  as  functions  of  airspeed  and  altitude,  however  these  mechaniza¬ 
tions  have  seme  shortcomings  which  are  pointed  out  in  this  paper,  and  also, 
the  desirability  of  self-adaptive  aircraft  control  systems  is  discussed. 

A  survey  of  adaptive  techniques  which  have  been  applied  to  the  airplane 
control  problem  is  presented,  and  a  self-adaptive  pitch-stability  aug- 
mentor  which  utilizes  a  digital  computer  in  the  adaptation  loop  is  pro¬ 
posed,  The  proposed  self-adaptive  system  is  tested  in  conjunction  with 
a  digital  three  degree  of  freedom  flight  simulation  of  the  pitch  dynamics 
of  a  representative  supersonic  fighter-attack  aircraft,  and  is  found  to  be 
capable  of  performing  the  adaptation  process  in  two  cycles  of  airplane 
short-period  oscillation.  Requirements  for  practical  mechanization  in 
a  flight  vehicle  and  recommendations  for  future  research  are  discussed. 
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As  the  speed  and  c  titude  capabilities  of  modern  high  performance- 
supersonic  aircraft  have  increased  in  recent  years,  the  range  of  varia¬ 
tion  of  the  dynamic  characteristics  of  these  aircraft  have  increased  pro¬ 
portionately.  A  ten  to  one  variation  in  damping  ratio  and  natural  fre¬ 
quency  of  some  of  the  dominant  modes  of  aircraft  motion  is  not  uncommon 
for  modern  airplanes  with  Mach  2  and  70,000  ft.  altitude  capabilities. 

As  an  example,  for  the  airplane  considered  in  this  study,  the  damping 
ratio  of  the  dominant  longitudinal  mode  varies  from  £  =  0.445  at  M  =  0.9 
@  10,000  ft.  to£  =  0.064  at  M  =  1.2  <?  45,000  ft. 

Various  stability  augmentation  devices  have  been  used  to  decrease  the 
range  of  variation  of  airplane  parameters  and  to  hold  them  as  near  as 
possible  to  optimum  values.  These  devices,  in  their  historical  order  of 
application,  include  fixed  gain  dampers,  pilot-adjusted-gain  dampers, 
speed-and-altitude- programmed -gain  dampers,  and  se If -adapt ive- gain -adjust¬ 
ing  dampers.  A  number  of  methods  have  been  investigated  for  mechanizing 
self-adaptive  stability  augmentation  systems  for  aircraft,  utilizing  vari¬ 
ous  techniques  which  will  be  discussed  later.  And,  the  purpose  of  this 
study  is  to  investigate  the  use  of  an  airborne  stored  program  digital 
computer  for  generation  of  the  adaptation  commands  in  a  self-adaptive 
pitch  damping  system  to  be  used  in  a  high  performance  fighter- interceptor 
type  aircraft.  The  dynamics  of  the  airplane,  including  control  system 
components,  and  the  adaptation  computer  were  simulated  on  a  Control  Data 
Corporation  1604  general  purpose  digital  computer. 

Before  beginning  the  discussion  of  adaptive  control  systems  as  applied 
to  aircraft,  it  is  felt  that  a  bettor  understanding  of  the  system  require¬ 
ments  may  be  gained  by  a  short  review  of  the  characteristic  dynamic  modes 
of  a  supersonic  aircraft  and  the  theory  underlying  the  analysis  and  syn¬ 
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The  airplane  ikj  i-j  axis  system  which  is  generally  used  in  airplane 
stability  .raid  c outre 3  studies  is  shown  in  figure  1,  as  three  mutually 
perpendicular  axes  intersecting  at  the  airplane  center  of  gravity.  The 
longitudinal,  or  X,  axis  lies  in  the  airplane's  vertical  plane  of  sym¬ 
metry  and  is  arbitrarily  oriented  by  the  airplane  designer  so  as  to  de¬ 
fine  a  horizontal  plane  of  symmetry  of  the  fuselage.  In  almost  all  air¬ 
craft  the  X  body  axis  does  not  coincide  with  the  longitudinal  principal 
axis  (longitudinal  axis  about  which  the  moment  of  inertia  is  least)  due 
to  the  high  location  of  the  mass  of  the  tail  surfaces  and  due  to  the 
assymetric  top-to-bottom  configuration  of  the  fuselage.  Consequently,  the 
longitudinal  principal  axis  normally  slopes  downward  from  tail  to  nose  and 
coincides  with  the  X  body  axis  only  at  the  center  of  gravity.  The  dynamic 
unbalance  of  mass  about  the  X  axis  thus  causes  the  airplane  to  have  an 
I  product  of  inertia  which  will  be  important  when  cross  coupling  of  air- 
plane  motions  is  discussed. 

The  Z  body  axis,  or  normal  axis,  is  perpendicular  to  the  X  axis  and  to 
the  wing  plane  and  has  positive  direction  downward.  And,  the  Y,  or  lateral 
axis,  is  simply  defined  as  being  normal  to  the  X  and  Z  axes.  Fortunately, 
due  to  the  lateral  symmetry,  an  airplane  has  no  dynamic  unbalance  about  the 
Y  axis  and  thus  only  one  product  of  inertia,  I  ,  exists. 

A  second  set  of  references,  the  earth  reference  system,  must  also  be 
considered  due  to  action  of  components  of  gravity  along  the  three  axes. 

To  determine  the  gravity  components,  the  airplane's  attitudes  in  pitch  and 
roll  must  be  determined  and  the  acceleration  of  gravity  resolved  along  the 
three  axes  by  utilizing  the  Euler  equations.  The  pitch  and  roll  angles  are 
illustrated  In  figure  2,  where  pitch  angle  is  show  with  wings  level  and 
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AXIS 

X 

y 

Z 

Units 

1.  Forces  acting  on  the  airplane 
along  the  axes: 

X 

Y 

z 

lb. 

2.  Moments  acting  on  the  airplane 
about  the  axes: 

L 

M 

N 

ft- lb. 

3. (  Moments  of  Inertia  of. the 
airplane  about  the  axes: 

lx 

iy 

'•Iz 

*  \ 

slug-ft^ 

4.  Linear  velocities  of  the 
airplane's  CG  in  space: 

u 

V 

w 

ft/sec 

5.  Angular  displacements  of  the 
airplane  about  the  axes: 

f" 

e 

rad. 

6.  Angular  velocities  of  the 
airplane  about  the  axes:.. 

p 

q  . 

.  J 

r 

rad/sec 

Figure  1.  Body  axis  system  and  parameters  used  in  the 
digital  simulation 
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reference  angles  a re.  equivalent  to  Integrations  of  airplane  angular  vel¬ 
ocities  about  their  respective  body  axes.  But,  if  the  craft  is  in  an 
attitude  with  non-zero  values  of  both  pitch  and  roll,  these  angles  must  be 
obtained  from  the  projection  of  the  X  and  Y  axes  onto  two  vertical  planes 
containing  these  axes, 

A  third  earth  reference  angle  which  has  not  been  illustrated  and  which 
plays  no  part  in  the  dynamics  of  motion  is  the  yaw,  or  heading,  angle  which 
is  the  projection  of  the  longitudinal  body  axis  onto  the  horizontal  plane. 

Characteristic  Airplane  Dynamics 

The  six  degrees  of  freedom  of  motion  of  an  aircraft  can  now  be  defined 
as  three  angular  freedoms  in  roll  pitch  and  yaw  about  the  X,Y*  and  Z  axes 
respectively;  and  longitudinal,  lateral,  and  normal  translations  along  the 
X,  Y,  and  Z  axes  respectively.  Thus,  the  dynamics  can  be  represented  in 
terms  of  six  simultaneous  differential  equations  of  motion.  Before  dis¬ 
cussing  the  makeup  of  the  six  equations  of  motion,  it  is  felt  that  a  better 
understanding  of  the  following  material  may  be  obtained  by  a  review  of  the 
various  modes  of  aircraft  motion.  The  motions  can  generally  he  divided  in¬ 
to  two  classes,  longitudinal  modes  and  lateral-directional  modes.  The  first 
class  involves  pitch  rotation  combined  with  translations  along  the  X  and  Z 
axes,  whereas  the  lateral-directional  modes  involve  rolling*  yawing^  and 
translation  along  the  Y  axis.  This  division  of  modes  occurs  because  of  the 
low  level  of  aerodynamic  and  inertial  cross  coupling  between  longitudinal 
and  lateral  motions.  For  small  amplitude  oscillations  (excursions  in  pitch, 
roll,  and  yaw  less  than  15-20  degrees)  the  two  classes  can  be  assumed  dynami¬ 
cally  independent . 
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Figure  2.  Earth  reference  system  used  in  the  digital  simulation 
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6  irscCior.nl  mo  las  of  a  rv  resentutlve  supersonic  a  Lrp  lane  <>  The  top  half 
of  the  S  plane  is  shown  for  each  case.  It  is  seen  that  the  longitudinal 
motions  occur  in  two  oscillatory  modes*  one  with  a  relatively  long  period, 
the  phngoid  mode,  and  one  with  a  much  shorter  period,  the  short-period 
mode.  The  damping  of  the  phugoic?  is  usually  lot*  and  in  some  cases  it  is 
a  neutral  or  even  divergent  oscillation.  However,  due  to  its  long  period, 
the  motion  is  easily  kept  under  pilot  control  with  little  or  no  ccxiscious 
attempt  at  correction,  and  no  artificial  stability  augmentation  is  re¬ 
quired  for  this  mode.  The  phugold  occurs  at  constant  angle  of  attack  and 
manifests  itself  as  a  continuous  trade-off  between  airspeed  and  alt  tude. 

Figures  4.  and  5.  show  that  the  period  and  damping  ratio  of  the  short 
period  mode  o£  a  supersonic  aircraft  can  vary  over  wide  limits,  and  the  in¬ 
clusion  of  the  time  lag  of  the  human  pilot  in  the  control  loop  can  some¬ 
times  cause  the  short-period  poles  to  move  into  the  right  hand  side  of  the 
S  plane.  This  necessitates  the  use  of  pitch  stability  augmentation  not 
only  for  improvement  of  handling  qualities  but  for  reasons  of  safety. 

This  mode  occurs  at  constant  airspeed  and  appears  as  a  pitching  oscilla¬ 
tion  accompanied  by  small  amplitude  translations  up  and  down  the  Z  axis. 

It  is  extremely  important  in  maneuvering  flight  and  is  of  prime  interest  to 
designers  of  airborne  fire  control  and  stability  augmentation  systems. 

The  lateral-directional  pole* zero  plot  for  a  supersonic  aircraft 
normally  exhibits  only  one  pair  of  complex  poles  and  tvc  real  poles  as 
shown  in  fig.  3.  The  frequency  and  damping  ratio  of  the  lateral-direction¬ 
al  short  period,  or  Dutch  roll,  mode  vary  as  flight  conditions  change  in  a 
manner  similar  tc  the  variation  of  longitudinal  short-period  characteristics. 
This  mode  is  a  combined  rolling-yawing  oscillation  with  small  amplitude 
lateral  translations,  and  its  name  is  derived  from  the  swaying  motions  of 
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Figure  3  Pole-Zero  Plot  of  the  Characteristic  Modes  of 
Motion  of  a  Supersonic  Aircraft. 
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a  Dutch  skater.  The  lutch  roll  oscillations  of  a  supersonic  aircraft 
usually  require  artificial  clamping  in  the  same  manner  as  required  by  the 
longitudinal  short-period  mode.  Therefore,  even  though  this  study  is  con¬ 
cerned  with  an  adaptive  longitudinal  stability  augmentation  system,  the 
results  and  conclusions  are  applicable  to  the  lateral-directional  case 
as  well. 

The  roll  damping,  or  roll  subsidence,  mode  is  shown  as  a  negative 
real  pole  and  it  represents  the  viscous  roll  damping  of  the  aircraft. 

And,  the  spiral  mode  is  usually  a  steady  divergence  which  manifests  itself 
as  a  gradually  steepening  turn  with  accompanying  sideslip.  The  divergence 
is  fortunately  slew  and,  as  in  the  case  of  the  longitudinal  phugoid  mode, 
is  easily  controlled  by  the  pilot.  The  roll  subsidence  and  spiral  mode 
poles  can  occasionally  form  a  complex  pole  pair  in  the  case  of  some  sub¬ 
sonic  airplanes  with  large  dihedral  effect.  This  lateral-directional 
"phugoid”  mode  will  exhibit  small  amplitude,  long  period  rolling  oscillations 
with  considerable  lateral  translation  and  very  little  yawing  oscillation. 

This  concludes  the  survey  of  normal  modes  of  airplane  motion.  How¬ 
ever,  there  are  other  "pathalogical"  modes  which  can  be  experienced  when 
the  rates  and  amplitudes  of  motion  become  large.  Some  of  these  are  the 
phenomena  of  inertial  cross  coupling  in  roll,  autorotation  (Ref.  125.)* 
and  longitudinal  pitch-up  (Ref.  86);  all  of  which  can  be  eliminated  or  al¬ 
leviated  by  the  use  of  well  designed  airplane  stability  augmentation  systems. 
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Aerodynamic  Stability  Derivatives 
One  final  topic,  the  concept  of  aerodynamic  stability  derivatives, 
should  be  mentioned  before  proceeding  to  the  equations  of  motion.  The 
forces  and  moments  which  act  on  an  aircraft  are  generally  expressed  in 
terms  of  non-dimensional  coefficients  multiplied  hv  appropriate  aircraft 
and  flight  condition  parameters.  As  an  example,  the  equation  for  airplane 
lift  is  shown  below, 

P  2 

Lift  *  CL  j  VS 

Where : 

Lift  is  in  pounds 

C  is  the  airplane  lift  coefficient 

Li 

p  s  atomspheric  density  in  slugs/ft~ 

V  55  airplane  forward  velocity  in  ft/sece 

2 

S  =  airplane  wing  area  in  ft  . 

P  2 

The  expression  —  V  is  a  measure  of  the  kinetic  energy  of  the  air- 
stream  and  is  termed  the  dynamic  pressure  because  it  also  represents  the 
force  which  would  theoretically  be  exerted  against  a  unit  square  surface 
held  normal  to  the  stream.  Thus,  the  product  of  dynamic  pressure  and  wing 
area  represents  a  hypothetical  aerodynamic  force  from  which  the  actual  air¬ 
plane  lift  can  be  obtained  through  multiplication  by  the  lift  coefficient. 
Moments  acting  about  the  aircraft  can  be  obtained  in  the  same  manner  with 
the  inclusion  of  a  moment  arm,  for  dimension  considerations,  as  is  done  in 
the  following  pitching  moment  equation. 

P  2  — 

Pitching  Moment  »  C  t  V  Sc 
°  m2 

Where  the  moment  arm  c  is  the  wing  mean  aerodynamic  chord.  Corres¬ 
ponding  expressions  for  drag  force,  side  force,  rolling  moment,  and  yaw¬ 
ing  moment  coefficients  are  included  in  the  table  of  symbols.  It  is  seen 
that  there  exists  an  aerodynamic  coefficient  for  each  of  the  six  degrees 
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of  freedom  or  motion,  These  coefficients  and  their  derivatives  are 
usually  determined  empirically  in  wind  tunnel  and  flight  studies  be¬ 
cause  of  the  inaccuracy  of  analytical  determinations. 

Each  coefficient  can  be  expressed  in  a  Taylor fs  series  such  as: 


Cr  m 

•••+(&) 


<Z+\dCC 


M  Ct.\  ao1 

\dCCLl  2.1 


■-H 


(c )CL\  . 

\7gTioc+  ladTaj 


_ oc  ‘ 

\dcc*J  2f 


+  ■ 


2! 


+ 


The  effects  of  the  second  order  and  higher  terms  are  insignificant  in  aero¬ 
dynamics  and  are  dropped.  So,  for  practical  purposes  the  coefficient  can 


be  expressed: 


c.-gD « *  6S 


OC 


-h 


(dCu 


The  partial  derivative  terms  such  as  are  called  aerodynamic  stability 

derivatives  and  are  usually  written  for  brevity  in  the  formC^^ 

•  .  CL-  CLccOO  d-  6lo>  CL  +  C iq 

Figure  6,  a  curve  of  lift  coefficient  vs,  angle  of  attack,  can  be  used  to 

illustrate  the  physical  meaning  of  stability  derivatives.  Here,(^-"r-) 

is  the  slope,  or  derivative,  of  the  airplane  lift  curve.  It  is  seen  that 

the  curve  is  linear  up  to  <X  ^  12°.  This  indicates  that  AC  ~  C CC  is 

L  hoc 

physically  valid  for  small  variations  in  OC  ,  The  fact  that  almost  all 
aerodynamic  stability  derivatives  remain  nearly  linear  over  a  similar  range 
of  their  corresponding  variables  enables  the  small  amplitude  motions  of  the 
aircraft  to  be  expressed  as  a  set  of  simultaneous  linear  differential  equa¬ 
tions. 

How  many  stability  derivatives  must  be  considered  in  the  equations  of 

■Ar  *«V 

motion9  Given  the  six  force  and  moment  coefficients  ,  Cxr>  ,  Ch  , 


*The  lift  coefficient  and  the  drag  coefficient  act  perpendi¬ 
cular  and  parallel,  respectively,  to  the  relative  wind,  i>ut  it  can  be 
assumed  with  little  error  that  they  are  also  normal  and  parallel  to  the 

X  body  axis.  Then  C_  =  C„  and  C  =  €_« 

Z  u  a  D 
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Angle  of  Attack,  deg. 

Figure  6.  Plot  of  Airplane  Lift  Coefficient  vs.  Angle  of  Attack. 

I 

and  C  and  the  fact  that  each  coefficient  can  be "inf luenced  by: 


1.  angular  displacement  about  each  axis, 

2.  angular  rate  about  each  axis, 

3.  velocity  along  each  axis,  and 

4.  acceleration  along  each  axis 


there  are  12  possible  derivatives  for  each  coefficient  or  72  in  all.  Also, 


control  surface  deflections  can  add  36  more  terms.  ^'Fortunately,  airplane 
symmetry,  which  causes  a  lack  of  cross  coupling  between  longitudinal  and 
lateral  -  directional  modes,  greatly  reduces  the  number.  The  magnitudes  of 
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some  of  the  remaining  derivatives  are  so  small  as  to  make  them  insigni¬ 
ficant,  As  a  result  there  are  approximately  25  derivatives  considered 
important  enough  to  he  included  in  a  simulation  program.  These  are  sum¬ 
marized  in  Appendix  TIT,  A  list  of  stability  derivatives  for  a  representa¬ 
tive  supersonic  aircraft  is  tabulated  in  Appendix  I.  The  aircraft  is  the 
same  one  for  which  the  longitudinal  short  period  mode  characteristics  are 
illustrated  in  figures  4,  and  5, 

The  Airplane  Equations  of  Motion  and  the  Longitudinal  Short  Period  Transfer 
Function, 

The  six  equations  of  motion  corresponding  to  the  six  degrees  of  free¬ 
dom  of  an  aircraft  are  derived  in  appendix  IV  and  the  final  results  are 
presented  in  table  I,  It  should  he  noted  that  the  equations  represent  the 
balance  of  inertial  and  aerodynamic  forces  and  moments  along  and  about  the 
airplane  body  axes,  and  that  the  terms  are  composed  of  .airplane  physical 
constants,  aerodynamic  stability  derivatives,  and  aerodynamic  variables. 
Assumptions  used  in  the  derivation  of  the  equations  are: 

1.  The  change  in  aerodynamic  forces  and  moments  due  to  deviations 
from  an  equilibrium  condition  are  assumed  directly  proportion¬ 
al  to  the  deviation. 

2.  The  deviations  of  the  airplane  from  equilibrium  are  assumed 
small, 

3.  The  airplane  is  a  rigid  body. 

4.  Constant  thrust  is  assumed  during  an  airplane  distrubance, 
and  secondary  power  effects  are  neglected. 

The  equations  are  used  directly  in  the  digital  computer  simulation  of  the 
airplane  dynamics,  which  is  outlined  in  appendix  VIS  and  the  airplane 
longitudinal  short-period  transfer  function  is  derived  from  the  pitching 
and  normal  translation  equations  of  motion  in  appendix  V.  This  transfer 


function  is  of  prime  importance  in  the  study  of  the  airplane  pitch  damper 
because  the  short -period  mode  is  the  dominant  longitudinal  mode.  The 
short-period  transfer  function  is  written  below  both  in  its  general  form  and 
also  in  terms  of  the  airplane  physical  characteristics,  airplane  stability 
derivatives,  ar.d  aerodynamic  parameters.  The  second  form  of  the  transfer 
function  is  shown  to  illustrate  the  dependence  of  the  dynamics  of  the  short 
period  mode  on  changes  in  stability  derivatives  and  aerodynamic  parameters, 
and  it  is  possible  from  inspection  of  the  function  to  determine  which  de¬ 
rivatives  and  parameters  exert  influences  on  a  given  characteristic,  such 
as  damping  ratio,  of  the  short-period  mode. 


Longitudinal  Short  Period  Transfer 


Function 


A.  General  Form 
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Ill  -  THE  !ttB  OF  A  FITCH  l  ':0  L0L40 IT UT) I hVJ.  ST.  mm  V  S  •LifTATIOH 

The  purpose  of  airplane  stability  augmentation  is  to  reluce  the 
range  of  variation  of  a  dynamic  characteristic,  or  characteristics*  and 
to  hold  the  variables  of  interest  ns  close  as  possible  to  predetermined 
optimum  values*  In  the  case  of  an  aircraft  these  optimum  values  are  de¬ 
fined  in  terms  of  "airplane  handling  qualities'*  which  describe  an  air¬ 
plane*  s  response  to  manual  control  by  a  pilot.  Even  though  the  determina¬ 
tion  of  desirable  handling  qualities  is  subjective,  and  is  dependent  on 
pilot  opinion,  a  number  of  studies  have  been  conducted  to  specify  these 
qualities  in  quantitative  engineering  terms*  One  of  the  better  known 
studies,  conducted  at  the  Cornell  Aeronautical  Laboratory,  (Ref*  34),  con¬ 
cludes  that  the  desirable  natural  frequency  and  damping  ratio  of  the  short- 
period  mode  of  a  fighter  airplane  lie  in  the  range  COn  s  .4  to  *5  cycles  per 
second  and  damping  ratio  ~  .65  to  ,70.  The  same  study  also  concludes  that 
even  though  the  longitudinal  natural  frequency  is  of  importance  in  maneuver¬ 
ing  flight,  it  can  vary  over  wide  limits  and  the  handling  qualities  will 
still  be  considered  tolerable,  as  long  as  the  short  period  damping  ratio  is 
held  in  the  .65  to  .70  range*  Thus  it  is  considered  more  desirable  to  equip 
a  high  performance  airplane  with  artificial  damping  than  to  provide  static 
stability  augmentation.  It  should  be  noted  that  for  the  craft  considered 
in  this  study,  as  well  as  for  almost  all  other  airplanes  of  similar  type 
and  mission,  the  longitudinal  damping  ratio  of  the  un-augmented  airplane 
is  always  below  the  desired  .7,  as  is  shown  in  figure  4.  This  means 
that  the  pitch  damping  system  will  work  to  supplement  the  basic  aerodynamic 
clamping  in  all  cases. 


18 


Mechanization  of  the  I  itch  hamper 


In  the  discussion  of  the  airplane  longitudinal  short “period  transfer 
function  it  was  shown  that  the  short-period  motion  can  be  represented  as  a 
second  order  system  with  characteristic  damping  ratio  and  natural  frequency. 
It  was  also  shown  that  the  natural  frequency  and  damping  ratio  can  be  ex¬ 
pressed  as  functions  of  the  aerodynamic  stability  derivatives.  The  ex¬ 
pression  for  short  period  damping  ratio  being: 


a  r  -  r 

mV  ^Lgc  £V7V 


From  the  above  relationship  it  can  be  seen  that  the  most  effective  means 
of  varying  the  damping  ratio  is  to  artificially  vary  the  pitch  damping  de¬ 


resistance  to  pitching  velocity. 

The  pitch  damper  must  then  provide  an  aerodynamic  pitching  moment  which 
is  proportional  to  airplane  pitch  rate,  and  which  has  a  direction  in  opposi¬ 
tion  to  the  angular  rate.  Thus,  the  moment  produced  is  in  phase  with  the 
viscous  aerodynamic  damping  in  pitch  and  the  two  moments  combine  to  form 
the  total  pitch  damping  moment.  The  artificial  damping  moment  is  usually 
generated  by  deflecting  the  airplane0 s  horizontal  stabilizer  in  direct 
proportion  to  measured  pitch  rate. 

Figure  7.  is  a  block  diagram  of  the  pitch  dampen  u ru:A  In  this  study* 

The  airplane  transfer  function  shown  in  the  diagram  represents  the  longi¬ 
tudinal  short -period  mode  only  because,  as  will  be  demonstrated  later,  the 
phugoid  mode  and  the  lateral-directional  dynamics  are  of  negligible  impor¬ 
tance  in  longitudinal  stability  augmentation  considerations. 

The  stabilizer  actuator  in  the  diagram  is  an  irreversible  hydraulic  position 
servo,  and  the  parameters  in  the  transfer  function  represent  the  state~o£- 
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Irreversible  PyJraulic  Airplane  Pit  on 
(Stabilizer  Actuator  Analc  Trccnsfer 

(Position  Qervo  )  Function 
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O 

Figure  7«  Block  Diagram  of  the  Pitch  D&mpar 
the-art.  as  used  in  modern  aircraft.  It  is  generally  assumed  that  the 
dynamics  of  the  rate  gyre  have  negligible  effect  on  system  operation  be¬ 
cause  the  ratio  of  the  characteristic  frequency  of  state-of-the-art  rate 
gyros  to  the  short-period  mode  frequency  is  of  the  order  of  20  to  30 
(Ref.  43),  and  thus  the  gyro  transfer  function  is  assumed  to  he  unify. 

Table  II  lists  the  II  airplane  transfer  functions  use  I  in  the  study  alon^ 
with  corresponding  flight  conditions.  And,  figures  S  thru  18  are  root 
locus  plots  showing  the  behavior  of  the  system  with  varying  dt-mp^r  gain  at 
each  condition.  The  square  symbols  on  the  plots  indicate  the  gair  values 
at  which  the  airplane  short-period  damping  ratio  has  been  increased  to  0.7. 
It  is  seen  that  the  use  of  the  carper  has  %*2ty  little  influence  on  the 
short  period  frequency  when  operate  \  over  the  range  of  gains  required  to 
bring  £  up  to  0.7  at  each  flight  condition*  And,  it  is  also  seen  that 


the  actuator 


pole  moves  i r,  such  a  direction,  with  increasing  damper  gain, 
that  the  actuator  damping  ratio  regains  close  to  0, 7  over  the  danper  opera¬ 
ting  range*  The  actuator  rode  will  become  unstable  at  extremely  high  gains* 
but  operation  at  such  gains  is  of  no  interest  in  this  study* 

It  is  known  that  the  closing  of  a  feedback  loop  around  a  plant  having 
varying  characteristics  will  cause  the  plant  output  to  vary  over  a  smaller 
range,  as  plant  characteristics  vary,  than  the  range  of  variation  of  out¬ 
put  with  no  feedback*  From  this  it  would  be  expected  that  the  pitch  damper 
should  not  only  Increase  the  average  value  of  damping  ratio  over  all  flight 
conditions,  but  also  that  the  range  of  variation  of  damping  ratio,  as  condi¬ 
tions  vary,  should  decrease*  This  is  borne  out  in  Figure  19.  which  is  a 
plot  of  short-period  damping  ratio  vs,  speed  and  altitude  with  constant 
damper  gain  s  0.3.  The  chosen  value  of  damper  gain  is  approximately  the 
median  of  the  gains  required  for  0*7  damping  ratio  over  all  11  flight 
conditions.  In  Che  case  of  the  basic  airplane,  £  varies  from  0.445  at  M 
=  0,9  <?  10,000  ft.  to  £  =  0.064  at  M  =  1*2?  45,000  ft.  whereas  the  varia¬ 
tion  of  /  vrith  damper  gain  »  0.3  is  from  1,1  to  0.436  for  the  same  two 
flight  conditions*  Thus  the  range  of  variation  about  the  mean  has  been 
reduced  from  approximately  190%  to  approximately  95%. 

Even  though  the  use  of  the  constant  gain  damper  has  caused  a  consider¬ 
able  improvement  in  the  dynamics  of  the  airplane,  a  further  restriction  in 
the  range  of  variation  of  damping  ratio  is  usually  desired  for  fighter-inter™ 
ceptor  type  aircraft*  One  method  of  solving  this  problem,  which  is  in 
general  use  today,  is  the  programming  of  stability  augmentation  system 
gains,  as  functions  of  measured  flight  parameters,  in  such  a  manner  that 
the  programmed  gains  will  cause  near  optimum  system  response  over  all 
flight  conditions.  The  techniques  of  gain  scheduling  methods  range  from 
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variation  oi  Jasper  gain  nr  function  of  e.casu red  airspeed  to  systems 
which  program  as  functions  of  airspeed*  altitude,  fuel  load,  center  of 
gravity  position,  wing  flap  position,  landing  gear  position,  and  vari¬ 
ous  combinations  thereof.  One  fact  becomes  obvious  from  the  preceeding 
discussion;  the  dynamic  characteristics  of  the  aircraft  must  be  accur¬ 
ately  known  in  advance  at  all  flight  conditions  and  for  all  airplane 
configurations  if  the  programming  system  is  to  operate  successfully  * 

Such  a  detailed  determination  of  airplane  characteristics  requires  a 
large  amount  of  wind  tunnel  and  flight  testing  which  is  both  time  consuming 
and  costly*  And,  even  after  the  desired  information  has  been  determined 
it  is  usually  necessary  to  conduct  lengthy  flight  tests  with  the  program¬ 
ming  and  stability  augmentation  systems  installed  in  the  airplane  for  final 
"tailoring*'  of  the  aircra ft- system  combination*  In  addition,  accurate  in¬ 
flight  air  data  and  airplane  configuration  Information  must  be  continusouly 
fed  into  the  system  during  flight  to  provide  the  basis  for  gain  scheduling* 
The  instrumentation  and  transducers  required  are  generally  costly,  and  are 
usually  complex  enough  that  they  create  reliability  and  fail-safety  problems, 
A  possible  solution  to  these  control  problems  is  to  design  a  system 
which  is  intelligent  enough  to  determine  when  the  airplane °s  dynamic  per¬ 
formance  has  changed,  and  in  what  manner  it  has  changed,  and  to  use  this 
information  to  adjust  the  stability  augmentation  system  gains  to  obtain 
optimum  performance  over  all  flight  conditions*  This  is  the  concept  of 
the  self-adaptive  control  system. 


33 


IV 


SELF-ADAPTIVE  CONTROL  SYSTEMS 


The  concept  of  adaptive  control  has  been  applied  to  a  wide  range  of 
control  system  types  with  the  result  that  a  number  of  definitions  of  adap¬ 
tive  control  exist.  However,  the  most  often  used  definition*  and  also  the 
one  which  appears  to  best  fit  the  aircraft  control  problem  is  as  follows: 

An  adaptive  system  is  a  control  system  configuration  in  which 
the  measurement  of  process  dynamics  or  signal  characteristics 
is  utilized  to  automatically  adjust  the  controller  in  an  attempt 
to  achieve  optimum  operation  at  all  times. 

From  this  definition  it  can  be  seen  that  the  adaptive  process  must  per¬ 
form  two  functions,  the  first  of  which  is  the  ,f identification”  function, 
or  the  measurement  of  process  dynamics  cr  signal  characteristics  to  deter¬ 
mine  present  system  state  and  compare  it  with  the  desired  state.  And,  the 
second  is  the  "command”  function  which  is  automatic  changing  of  the  system 
so  that  it  will  approach  a  desired  state.  In  most  adaptive  control  studies 
it  is  assumed  that  the  variation  of  system  characteristics  takes  place  slov; 
ly  with  respect  to  the  frequency  of  the  system  mode  of  interest,  and  thus 
the  identification  process  can  be  allowed  a  relatively  long  measurement  in¬ 
terval  during  which  to  determine  system  state.  However,  in  the  ail  craft 
problem,  it  is  desirable  to  perform  the  adaptation  process  relatively  rapid 
ly  in  comparison  to  the  quantity  being  controlled,  in  this  case  the  short 
period  mode  dynamics.  This  is  because  virtually  instantaneous  changes  can 
take  place  in  the  airframe  transfer  function  due  to  release  of  payload  or 
certain  Mach  number  effects  on  stability  derivatives.  For  this  reason* 
one  of  the  objectives  of  this  study  is  the  investigation  of  an  adaptive 
technique  capable  of  rapid  identification  and  command.  Another  considera¬ 
tion  in  mechanization  of  the  adaptive  process  is  the  generally  accepted 


requirement  feu*3  >  iu  o  :Iei  :oi  the  system  it  be  con.si.deve  1  truly  adap- 
tive,  the  identif itaiion  process  can  net  involve  measurement  of  any 
variable  which  Is  external  to  the  system.  To  illustrate  the  concepts 
and  requirements  discussed  above,  a  block  diagram  of  Gibson* s  ideal 
generalised  adaptive  system  (Ref.  76)  is  shown  in  Figure  20 . 


Identification  and  Co  mound 
Computer 


CONTROLLER 


MEMORY 


COMPUTATION 


COMMAND 

GENERATION 


Comoran  d  SioncL !  for 


Control ier  /Alteration 


Figure  20,  Gibson1 s  Ideal  Generalized  vlaptive 

System 

Classification  of  Adaptive  Systems 
In  o'cdei'  to  compare  the  approaches  which  have  been  taken  to  the 
design  of  adaptive  control  systems,  Aseltlne,  in  reference  4*  has  seperat 
ed  them  into  5  classes  which,  in  some  cases,  have  a  certain  amount  of 
overlapping,  but  even  so,  they  form  a  convenient  framework  for  comparison 
Aseltine's  classes  are: 
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1)  assive  Vi  i  tat: ion:  Fastens  winch  adapt  without  automatic 
changes  ir»  system  characteristics,  but  rather  through  de~ 
sign  for  uniform  operation  over  wide  variations  in  environ¬ 
ment  . 

2)  Input  Signal  Adaptation:  Systems  which  automatically  ad¬ 
just  their  parameters  in  accordance  with  input  signal 
characteristics* 

3)  Sxtrenam  Adaptation;  Systems  which  adjust  to  obtain  a  maxi¬ 
mum  or  minimum  of  some  system  variable, 

4)  System  Variable  Adaptation:  Systems  in  which  the  command 
process  is  actuated  directly  by  measurement;  of  a  system 
variable  or  variables* 

5)  System  Characteristic  Adaptation:  Systems  in  which  the 
identification  process  measures  system  variables  in  order 
to  determine  system  characteristics.  The  system  character¬ 
istics  thus  determined  are  compared  with  desired  character¬ 
istics  and  the  differences  used  to  generate  command  signals* 

It  has  been  noted  from  a  survey  of  adaptive  systems  which  have  been  invest 
gated  for  possible  use  in  aircraft  control  that,  to  date,  all  such  systems 
fit  into  classes  1,  4,  or  5.  And,  it  has  also  been  noted  that  an  alter* 
nate  classification  system  might  be  applied  to  aircraft  adaptive  control 
systems.  The  suggested  classes  are: 

1)  High  Gain  Adaptation:  System s  in  which  the  characteristics 
of  the  aircraft  are  submerged  by  the  use  of  high  feedback 
and/or  forward  loop  gains.  Some  form  of  Idealized  model  is 
normally  used  to  provide  either  system  input  signals  or  out¬ 
put  comparison  signals,  and  the  techniques  usually  involve 
either  nonlinear  processes,  with  attendant  limit  cycle,  or 


operation  of  a  minor  system  mode  near  instability. 

2)  Impulse  Response  Adaptation;  Systems  in  which  the  identi¬ 
fication  process  computes  a  time  history  o£  the  plant  im¬ 
pulse  response.  The  impulse  response  is  then  utilized  by 
the  command  process  to  generate  plant  alteration  signals, 

3)  Parameter  Adaptation:  Systems  in  which  plant  parameters 
are  determined  by  direct  computation.  The  error  between 
the  desired  value  of  the  parameter  and  its  actual  value  is 
used  as  the  command  signal  for  plant  alteration. 

Survey  of  Adaptive  Techniques  which  have  Been  Applied  to  the  Airplane 

Stability  Augmentation  Problem 

1)  High  Gain  Adaptation  Systems 

The  earliest  studies  of  aircraft  adaptive  control  involved  systems 
which  operated  at  high  gains  in  order  to  force  the  variation  of  airframe 
parameters  into  a  small  range.  One  of  the  first  of  these  studies  was 
conducted  by  Rath  of  WADC  (Ref.  134)  on  a  method  which  was  not  truly  ad¬ 
aptive,  in  that  parameters  were  not  automatically  varied*  however*  it 
may  be  classed  as  passive  adaptation.  Rath* s,  system  used  a  parallel  com¬ 
plementary  optimum  response  model  as  shown  in  Figure  21.  The  goal  of 
this  study  was  to  design  a  system  which  would  produce  optimum  airframe 
dynamic  characteristics,  as  airframe  parameters  varied,  with  constant 
system  gain.  A  fourth  order  transfer  function  was  selected  as  the  opti¬ 
mum  model  for  the  pitch  dynamics  and  the  same  input  signal  was  applied  to 
the  airframe  and  to  the  model.  The  difference  between  model  and  plant 
output  was  obtained  and  various  functions  of  the  difference  signal  were 
utilized  to  improve  the  airframe  response,  Both  the  forward  loop  gain 
and  the  shaping  loop  gain  were  made  high  in  order  to  M force”  the  system 


Figure  21.  Rath*s  System  for  passive  Adaptation  Using  a  Complementary 

Optimum  Response  Model 

overall  characteristics  to  resemble  the  model.  Feedback  functions  studied 
included  pitch  angle*  or  position*  error  feedback,  integrated  position 
error  feedback*  rate  error  feedback*  and  rate-error-plus-position-error 
feedback.  It  should  be  noted  that  error*  as  referred  to  here*  means  the 
difference  between  plant  output  and  model  output  rather  than  plant  input- 
output  error.  The  first  two  feedback  functions  produced  overall  instabil¬ 
ity  at  high  gains*  whereas  rate  error  feedback  resulted  in  a  stable  system 
with  infinite  steady  state  position  error.  Rate-error-plus-position-error 
feedback  allowed  very  high  gains  to  be  used  before  the  system  became  un¬ 
stable*  and  while  the  system  was  stable*  the  airframe  output  rr&tcbed  the 
model  output  acceptably  at  most  flight  conditions  while  using  a  constant 
value  of  feedback  gain.  It  was  also  found  that*  whatever  feedback  shaping 
was  used,  the  basic  airframe  and  hydraulic  actuator  must  fom  a  stable 
system  with  the  shaped  feedback  loop  alone.  That  is*  the  feedforward  loop 


could  rot  coft  cc  m-st  *  >L  {■  >  of  rhe  a.-.’  shaped  fV^Jhack 

combination.  The  optimum  **ra  *?fer  function  feedforwrd  in  combination 
with  the  shape!  feedback  h**d  the.  effect  of  changing  the  input  to  the 
airfrai  a  to  a  value  which  caused  'the  airframe  to  have  the  desired  out¬ 
put,  assuming  the  airframe  htl  th*  desired  control  rates  and  high  gain 
control  power « 

An  approach  similar  ro  RairhAs  passive  adaptation  Method  was  taken  by 
Jandols  of  Corwair  who  utilised  the  inverse  of  the  optimum  response  model 
in  the  feedback  loop  in  combination  with  a  high  forward  loop  gain  system 
as  shown  in  Figure  22*  (Ref*  43). 


AIRFRAME  SHORT  PERIOD 
TRANSFER  FUNCTION 


Figure  22*  Dandois 8  System  Using  the  Inverse  of  the  Optimum  Response  Model 

in  the  Feedback  Loop 

He  demonstrated  that,  with  such  a  feedback  arrangement,  the  overall  system 
characteristics  approached  the  characteristics  of  the  optimum  response 
model  as  forward  loop  gain  was  increased*  In  actual  operation,  the 
system  held  airrrana  parameters  wichir  accet  table  lir.its  c»&r  -  wide 


range  oC  flight  con ilt ions  and  the  i\ghea.  the  gain  the  better  the 
performance,  However,  gain  1  units,  dependent  upon  "light  conditions, 
were  reached  at  which  the  systen  became  unstable.  The  instability  was 
shown  to  be  due  to  sensor  non linearities  rather  than  a  basic,  theoretical, 
system  instability. 

One  of  the  first  "high  gain"  adaptive  systems  to  actually  use  an 
identification  and  command  process  was  the  system-variable-adaptation 
system  of  Fluggc-Lotz  and  Taylor  (Ref,  64),  which  is  Illustrated  in  Fig* 
ure  23,  As  is  shown  in  the  block  diagram,  one  of  a  number  of  alternate 


Figure  23,  Nonlinear  Adaptive  System  of  Flugge-Lotz  and  Taylor 
feedback  paths  was  used  according  to  a  switching  logic  process  based  on 
simultaneous  measurements  of  system  output,  error  signal,  and  derivatives 
of  these.  The  Identification  process  thus  consisted  of  the  measurements 
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*n>  sv  ■  tch*  "a  : «.  ••*»  u  •  rb  •  *.pj.  rvf  u/tb  luaU  path,  and 

the  coniP-.ati^  cons  ■  sle  of  in-,  itching  to  the  chosen  path.  Forward  gain  was 
kept  high  in  order  to  “force”  a  small  instantaneous  error  between  input 
and  output,  and  the  system  was  nonlinear  during  any  length  of  time  com-* 
parable  to  the  natural  period  of  the  plant  because  of  the  switching  be¬ 
tween  feedback  loops  during  one  cycle  of  operation*  The  arrangement  was 
shown  to  be  insensitive  to  variations  in  plant  parameters  over  wide  limits, 
however  a  limit  cycle  persisted  during  operation  due  to  the  feedback  non¬ 
linearity. 

The  Minneapclis-Honeyvel!  relay  actuated  adaptive  system  (Ref*  165), 
which  is  presently  undergoing  flight  tests,  is  similar  to  the  Flu&ge-hotz 
mechanization  but  uses  only  one  feedback  path  with  switching  logic  in  the 
forward  loop.  The  switching  logic  derives  Its  input  from  the  error  signal 
between  optimum  response  model  output  and  the  process  feedback  signal  as 
shown  in  Figure  24*  The  switching  logic  actuates  a  relay,  or  on-off  con 
troller,  to  operate  the  process  at  a  high  gain  with  the  objective  of 
maintaining  a  small  error  signal  between  model  output  and  plant  output* 

As  in  the  Flugge-Lots  system,  the  airframe  output  approximately  dupli¬ 
cates  the  model  output  when  the  forward  gain  is  high,  and  the  character¬ 
istic  limit  cycle  is  also  present.  A  second  adaptive  feature  of  this 
system,  aside  from  the  switching  logic,  is  the  us^  of  limit  cycle  ampli¬ 
tude  measurement  to  vary  forward  loop  gain.  That  is,  forward  gain  is  ad¬ 
justed  in  such  a  manner  that  the  limit  cycle  amplitude  is  held  constant 
over  all  flight  conditions.  The  prototype  of  this  system  was  successfully 
tested  in  an  F-94C  fighter  aircraft  and  an  improved  version  is  presently 
being  flight  tested  in  the  X-15  high  altitude  research  airplane. 

Another  application  of  the  utilization  of  the  inverse  optimum  re¬ 
sponse  transfer  function  ixi  the  feedback  loop,  in  conjunction  with 


Limit  Cyo/e 


Amplitude  Meats  ur&noeni 


ACTUATOR 


|  RATE 

\GYRO 

Figure  24.  MJnneapolis-Moneywe! 1  Relay  Actuated  Adaptive  System 
high  forward  loop  gain,  is  the  General  Electric  system  (Ref.  32),  which 
was  recently  tested  in  an  F-102  interceptor  aircraft.  The  mechanization 
is  similar  to  that  of  Dandois  shovn  in  Figure  22,  except  that  ar>  active 
system-variable-adaptat ion  feature  has  been  added.  This  feature  Involves 
variation  of  the  forward  loop  gain,  as  flight  conditions  change*  to 
maintain  the  system  just  below  the  boundary  of  instability.  As  was  pro- 
viously  mentioned,  when  considering  systems  of  this  type*  the  higher  the 
gain,  the  closer  will  be  the  airframe  characteristics  to  the  optimum  re¬ 
sponse  characteristics  and  thus  if  the  system  Is  operated  near  the  in¬ 
stability  boundary, the  maximum  performance  will  he  extracted.  In  opera¬ 
tion,  the  adaptation  ^  stem  monitors  the  frequency  the  control  surface 
actuator  because  the  actuator  mode  is  the  one  which  becomes  unstable  as 
forward  gain  is  increased.  The  g*in  Is  then  adj  isted  to  iraintn  Ic  the 
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of  this  rode  ecus 


It  should 


frequency  of  this  mode  constant  over  a  1 1  flight  conditions, 
be  noted  that:  the  actuator  mode  is  a  lightly  damped  oscillation  and 
that  normal  control  Inputs  are  utilized  to  provide  the  required  pertur¬ 
bations.  The  frequency  measurement  is  accomplished  with  a  narrow  band 
pass  filter  and  zero  crossing  counter  arrangement. 


2}  Impulse  Response  Adaptation  Systems 

The  function  of  the  identification  process  In  the  impulse  response 
adaptation  class  of  systems  is  first,  the  determination  of  a  time  history 
of  the  plant  impulse  response  through  measurement  of  system  input  and  out¬ 
put  signals  over  a  period  of  time,  and  second,  the  determination  of  plant 
characteristics  through  examination  of  the  impulse  response.  An  early 
method  of  impulse  response  self  optimization  which  was  investigated  by 
Aseltine  (Ref.  4}  serves  to  illustrate  the  process,  even  though  this 
particular  mechanization  would  not  be  suitable  for  aeronautical  use. 

The  impulse  response  was  obtained  by  periodic  input  pulses  which  were 
widely  spaced  with  respect  to  the  period  of  the  plant  mode  of  interest. 
After  each  perturbation,  the  system  damping  ratio  was  computed  by  separ¬ 
ate  integration  of  the  positive  and  negative  portions  of  the  inpulse  re¬ 
sponse  and  obtaining  the  ratio  of  positive  to  negative  area.  Damping 
ratio  was  computed  from  ibis  ratio  and  the  system  damping  coefficient  was 
varied  in  a  closed  loop  manner  to  cause  actual  damping  ratio  to  approach 
a  desired  value.  This  met he  1  was  considered  undesirable  for  aircraft  use 
primarily  because  normal  operating  signals  must  be  excluded  from  the  in¬ 
put  during  the  comparatively  lengthy  identif icatior  period. 

A  method  of  impulse  response  determination  in  which  normal  Input 
signals  could  be  included  In  the  identification  process  was  investigated 


by  / nee r son  c  «:  h.ivj?  -in’  s  ^  ’  »  fho  *ech*vii.?ai;  ion  «ane  use  of  the 

convolution  theorem  '.-heuzby  the  system  input  convolved  In  t ini';  with  the 
impulse  response  lesuHs  in  a  time  history  of  the  system  output*  In  opera¬ 
tion,  random  noise  'with  auto*:  orrel  sti  cm  Z*mrtto\\(pii('Y)  was  superimposed  on 
the  norma!  input  signal  is  shown  in  Figure  25.  The  plant  output  was  filt¬ 
ered  to  separatee  the  noise  component  front*  the  tot  a  J  resj  o>nse»  and  the  in- 
put  and  output  noise  ccrx.t.  uents  were  cross-correlated  to  obtain  tf)jo (T)  - 

CO 

.  The  noise  input  had  a  bandwidth  3  to  10  times  the  ph^nt 

co 

bandwidth  so  that  ^ f)ii(T ')  was  af  fectively  impulse  and  (pio(/)  an  - 

proximate ly  the  plant  impulse  response.  Twelve  channels  of  corss  cor re la- 
lion  were  used  and  e^ch  channel  provided  one  point  on  a  time  history  of 
the  impulse  response.  Damping  ratio  was  determined  an:I  utilized  as  its  - 
cussed  in  the  proceeding  paragraph. 


PARAMETER  PARAMETER 


Figure  25. 


Block  Diagram  cf  :ercontronic*s  Impulse  Ttesporse  ce]f  {Vtimi-  - 

i vg  System 


/  e 


Another  mo  the  »  oc  b.«t  1  se  res;  onr>t  let*  rminauirm  the  convolu¬ 

tion  theorem  in  conjunction  with  simultaneous  input  and  output  measurement s 
is  currently  under  investigation  at  Hughes  Aircraft e  In  this  method,  d i s - 
crete  sampler  the  input  And  output  time  histories  are  use!  as  coeffici¬ 
ents  in  two  polynomials.  The  output  polynomial  is  periodically  divided  by 
the  input  polynomial  through  the  use  of  a  high  speed  digital  computer.  The 
coefficients  of  the  resultant  polynomial  can  be  shown  to  be  amplitude  sam¬ 
ples  of  the  system  impulse  response.  One  requirement  for  this  method  is 
that  the  amplitudes  of  the  system  output  samples  must  be  corrected  for 
initial  conditions  existing  at  the  beginning  of  the  interval  over  which 
samples  are  taken.  The  correction  process  requires  a  knowledge  of  the 
impulse  response  and  thus  the  impulse  response  must  be  determined  in  a  con¬ 
vergent  succession  of  polynomial  divisions.  At  present,  the  requirements 
for  convergence  to  the  true  impulse  response  are  under  investigation. 

3)  Parameter  Adaptation  Systems 

The  function  of  the  identification  process  in  n  parameter  adaptation 
system  is  the  direct  determination  of  plant  parameters  through  measurements 
of,  and  operations  upon,  selected  signals  and  their  derivatives. 

One  of  the  first  studies  of  parameter  adaptation  was  conducted  by  Kal¬ 
man,  of  Columbia  University,  who  used  a  digital  computer  to  determine  plant 
parameters  (Ref.  91).  In  Kalman8 s  system,  coefficients  of  the  plant 8 s  pulse 
transfer  function  were  computed  by  continuously  measuring  and  storing 
sampled  values  of  plant  input  and  output.  After  a  series  of  samples  had 
been  stored,  values  of  the  transfer  function  coefficients  were  assumed  and 
a  series  of  sample  output  points  were  digitally  calculated  using  the  known 
past  input  samples.  The  mean  square  error  between  calculated  and  known  out¬ 
put  points  was  determined  and  the  assumed  values  of  the  coefficients  were 
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alter-  *  hrc-ryt  r}  o  >•  1«ast  squ*  r*  s  which 

changed  coefficients  in  n  manner  that  reduced  the  near*  square  error.  This 
process  was  repeat'.*-1  At  very  high  s\  eed  until  the  MTT  was  reduced  be  lev  a 
desired  limit.  After  the  plant  coefficients  had  been  determined  with  suf¬ 
ficient  accuracy,  an  optimal  controller  was  synthesized  in  a  second  digital 
computer  routine  and  the  results  were  used  to  adjust  the  system  controller. 

A  different  approach  to  parameter  adaptation  was  taken  by  Whitaker,  of 
MIT,  who  used  a  model  reference  as  shown  in  Figure  26.  (Ref.  171). 


MEASURED  OUTPUT  QUANTITIES 
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Figure  26.  MIT  Model  Reference  Parameter  Adaptation  System 
The  Input  signal  was  fed  to  both  the  plant  and  to  the  ideal  model;,  and  a 
difference  function  between  plant  output  and  model  output  was  used  to  form 
an  integral  squared  error.  The  ISE  thus  formed  was  required  to  be  related 
to  a  single  plant  parameter,  which  was  then  varied  in  a  closed  loop  manner 
to  cause  the  ISE  to  attain  its  minimum  value.  The  primary  objective  of  the 
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study  was  the  determination  of  what  functions  of  the  plant-output- to-model- 


output  error  were  indicative  of  what  parameter  variations.  After  the  de¬ 
sired  relationships  were  determined,  the  system  was  mechanized  and  flight 
tested  in  an  F-94A  fighter  aircraft  and  it  was  found  that  adpatation  could 
occur  within  an  interval  corresponding  to  two  or  three  time  constants  of  the 
dominant  system  response  modes.  However,  recent  developments  have  made 
possible  a  significant  reduction  in  convergence  time  and  the  elimination 
of  much  of  the  complexity  of  the  previous  mechanization  (Ref.  127).  These 
developments  were  derived  from  a  recent  doctoral  research  y  rograrn  at  MI? 
which  evolved  an  analytical  design  method  for  parameter  adaptation  model 
reference  systems  * 

A,  method  of  adaptation  employing  an  analog  computer,  which  calculated 
values  of  instantaneous  forward  loop  gain,  was  studied  by  Corbin  of  Brook¬ 
lyn  polytechnic  Institute  (Rtf.  38).  The  gain  value  calculated  was  com¬ 
pared  with  a  desired  level  and  the  difference  used  as  a  signal  to  the  gain 
controller.  The  differential  equations  of  motion  of  the  plant  were  set  up 
in  the  analog  computer  as  integral  equations  and  were  solved  for  the  gain 
value,  given  continuous  measurements  of  the  plant  input  and  output  signals, 


as  shewn  In  Figure  27. 


Figure  27.  Corbin* r  Gain  Adaptation  Method  Utilizing  an  Analog  Computer 

in  the  Identification  Irocess 


A  method  similar  to  Corbin9 s  Approach,  but  which  employed  a  digital 
computer  to  directly  solve  for  the  plant  parameters,  was  investigated  by 
Azgapetian  of  Servomechanisms  Incorporated  (Ref,  7*)*  In  ibis  study,  the 
equations  of  motion  of  the  second  order  system  under  investigation  were 
converted  to  sample  data  notation,  and  relationships  for  damping  ratio, 
natural  frequency,  and  forward  loop  gain  were  derived  in  terms  of  past 
sample  points  of  plant  input  and  output.  The  equations  were  then  mechaniz¬ 
ed  in  a  digital  computer  routine  which  periodically  sampled  the  input  and 
output  signals  and  computed  plant  parameters  at  a  high  rate.  Parameters 
thus  calculated  were  used  in  a  closed  loop  manner  to  vary  a  controller  which 
caused  overall  system  dynamics  to  remain  constant  as  plant  characteristics 
were  varied. 

The  adaptation  method  investigated  in  this  study  is  similar  to  that 
discussed  in  the  preceeding  paragraph  in  that  a  digital  computer  is  util¬ 
ized  to  solve  second  order  equations  of  notion  for  system  damping  ratio 
and  natural  frequency  given  system  output  samples.  The  primary  differences 
are  the  method  of  solution,  and  the  fact  that  the  system  under  investigation 
is  actually  a  sixth  order  system  which  includes  a  second  order  hydraulic 
servo  and  a  fourth  order  airframe  transfer  function. 
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V  -  I^OFOSED  SYSTEM 

As  mentioned  in  the  preceding  section,  It  is  proposed  that  an  air¬ 
craft  pitch  damping  system  could  be  made  to  produce  invariant  short  period 
damping  ratio  over  all  flight  conditions  by  digital  computation  of  the 
short-period  damping  ratio  and  closed  loop  adjustment  of  the  pitch  damper 
gain  to  cause  the  damping  ratio  to  equal  a  desired  value,  A  block  diagram 
of  the  proposed  system  is  shown  in  Figure  28  below„  The  proposed  identifi¬ 
cation  process  thus  consists  of  solving  the  airframe  equations  of  motion 


Figure  28,  F reposed  Self-Adaptive  Fitch  Damper 
for  values  of  damping  ratio  given  samples  of  pitch  angle  and  its  deriva¬ 
tives,  and  the  command  process  consists  of  varying  damper  gain  at  a  rate 
proportional  to  the  error  between  desired  and  computed  /  and  in  the 
direction  indicated  by  the  sign  of  the  error,  Basic  assumptions  and  condi¬ 
tions  for  this  study  are: 
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1.  The  airframe  and  control  system  are  linear.  This  has 
been  shown  to  be  true  for  small  perturbations « 

2.  There  is  no  control  surface  rate  limiting  or  deflection 
limit ing. 

3.  Airframe  parameters  are  time  variant  but  are  assumed  con¬ 
stant  over  the  length  of  time  required  for  one  identifica¬ 
tion  cycle. 

4.  System  disturbances  due  to  external  noise  (atmospheric 
turbulence)  will  not  be  considered. 

It  is  also  assumed  that  the  short-period  mode  is  dominant  in  the  longitudin¬ 
al  dynamics  and  thus  the  airframe  longitudinal  response  can  be  closely  des¬ 
cribed  by  the  short-period  frequency  and  damping  ratio.  This  assumption 
was  verified  by  conducting  two  longitudinal  flight  simulation  programs  at 
each  of  the  11  conditions  listed  in  Table  II,  The  first  program  was  run 
in  two  degrees  of  freedom,  pitch  and  normal  translation,  and  thus  simulated 
the  short-period  mode  dynamics,  and  the  second  program  simulated  the  short- 
period  and  phugoid  modes  by  inclusion  of  the  longitudinal  translation  degree 
of  freedom.  Damping  ratio  and  period  were  determined  from  pitch  angle  time 
histories  obtained  by  pulsing  the  aircraft  about  the  pitch  axis  and  al loving 
it  to  return  to  equilibrium.  Characteristics  obtained  from  the  simulation 
programs  are  listed  in  Table  III,  It  is  seen  that  the  two  degree  and  three 
degree  of  freedom  characteristics  are  in  close  agreement,  and  thus  the  deci¬ 
sion  was  made  to  investigate  a  computational  method  of  determining  the 
'•pseudo"  second  order  damping  ratio  of  a  sixth  order,  three  degree  of  free¬ 
dom,  airframe-control  system  combination  from  sampled  data. 
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Table  III  -  Short -I  ax  10J  Ho  ’e  Characteristics  Cora;  are!  with  Dairying  uatio 
arc’  leriod  Obtained  froi..  a  longitudinal  Three  Degree  of  Ftee- 
dora  Flight  Simulation. 


Flight 

Condition 

Jet  in^ 

Rat  io 

Period, 

See 

M  '■ 

A1 t ItuCe 

2  DF 

3DF 

2DF 

3DF 

0.9 

lO.CCO1 

0.444 

0.440 

1.650 

1.700 

1, 2 

10, 000 1 

0. 135 

0.177 

1 .050 

1 . 005 

1.4 

10, 000 1 

0.155 

0. 156 

1.025 

1.020 

0.-9 

25,000 

0.245 

0.244 

1.974 

2.025 

1.2 

2.5,  roc1 

0.100 

C. 1023 

1.225 

1.275 

1.4 

2 5,0 oo 1 

0.104 

1.105 

1.199 

1 . 175 

1.3 

25,00c1 

0.095 

0.096 

1.200 

1.225 

2.9 

45, 000 1 

0. 140 

0. 141 

2.900 

2.075 

i  n 

t.  t  — 

45,0  00 1 

0.004 

0.066 

1.725 

1,725 

1.4 

45, 000 1 

0.063 

0.067 

1.650 

1.675 

1.3 

45,00c1 

0.057 

0.C58 

1 . 600 

1.599 
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The  generalized  differential  equation  of  motion  for  a  second  order 
system  is: 

Q  +  LSU)ne+W£G=  Ku 

from  which  an  expression  for  damping  ratio  can  be  derived  as. 

/=  ,  _  o  ~  6°^6 

7  zujne  2.Q 

It  is  seen  that,  in  this  form,  the  damping  ratio  is  a  function  of  the  in¬ 
stantaneous  values  of  8  and  its  first  and  second  derivatives,  the  input 
forcing  function,  and  the  natural  frequency.  Azgapetian  has  shown  in  Ref. 
7.  that  an  expression  for,  £  c<^ich  is  independent  of  CJn  ,  can  be  de¬ 
rived  as  an  integral  equation.  But  it  was  decided,  for  this  study,  to 
avoid  long  term  integrations  and  tueii  attendant  drifts,  and  rather  to  make 
use  of  the  high  speed  iterating  capabilities  of  digital  computers  by  de¬ 
riving  a  companion  expression  for  natural  frequency  to  he  used  in  an  itera¬ 
tive  process  for  arriving  at  £  and  U)n  simultaneously.  The  expression 
for  natural  frequency  being: 


U=  -  L§—  + 
n  e 


+ 


K§e 

e 


which  is  in  terms  of  instantaneous  values  of  0  and  its  first:  two  deriva¬ 
tives,  control  input,  and  damping  ratio.  It  is  noted  that  the  term,  k, 
by  which  the  control  deflection  is  multiplied,  corresponds  to  the  airplane 
elevator  effectiveness,  or  steady  state  airframe  gain,  which  is  also  a  func¬ 
tion  of  flight  conditions  and  which  therefore  must  be  determined  in  a 
spearate  computation  routine.  Corbin,  in  Ref.  38.  has  devised  a  method 
of  computing  the  instantaneous  value  of  forward  loop  gain  of  a  second  order 
system  by  utilizing  integral  equations.  Corbinfls  system  utilized  an  ana¬ 
log  computer  which  was  supplied  continuous  values  of  plant  input  and  out¬ 
put,  but,  the  method  could  be  converted  to  digital  form  with  sampled  in¬ 
puts  and  used  in  this  study  if  desired.  It  was  decided  not  to  introduce 
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this  additional  comnlexiM  irto  the  pro gran  but  rather  to  concentrate 
solely  on  the  plant  parameter  computation  method  and  its  application 
to  an  adaptive  system  because  the  computation  method  was,  as  yet,  untried. 
This  decision  necessitated  the  exclusion  of  pilot  control  inputs  to  the 
system  because,  if  pilot  inputs  were  to  be  included,  the  value  of  the  air¬ 
frame  gain  would  have  to  be  known.  It  should  be  made  clear  that  control 
deflections  called  for  by  the  pitch  damper  are  not  included  in  the  preceedin 
expressions  for  natural  frequency  and  damping  ratio  because  the  airplane 
pitching  moment  created  by  the  damping  system  is  considered  to  be  added  to 
the  natural  aerodynamic  damping  moment  to  make  up  total  damping.  If  pilot 
inputs  were  to  be  included,  the  control  deflection  created  by  the  pitch 
damper  would  be  subtracted  from  total  deflection  and  the  remaining  control 
deflection  would  be  used  as  the  input  term,  Sg  ,  in  the  /  and  CO n  equa¬ 
tions,  With  pilot  control  input  excluded,  the  equations  are  reduced  to: 

-^±  f 

e _ oo n  9 

ZLJnG  £0 

The  proposed  iterative  computational  routine  is  as  follows: 

1)  Measure  and  store  sampled  values  of  0,  9,  and  0  at  two 
separate  instants  in  time, 

2)  Assume  a  value  of  damping  ratio  and  compute  using  the 

first  sample  set  of  0,  0,  &  9  and  assumed  .  The  CJn 

equation  is  used  in  the  initial  calculation  because  the 
value  of  0)r)  is  dependent  upon  variations  in  £  to  a  lesser 
extent  than  is  the  value  of  £  dependent  upon  variations  in 
COn  . 

3)  Use  the  value  of  U)n  computed  in  2)  and  the  second  sample  set 
•  #  # 

of  9,  9,  &  ©  to  compute  £  . 


4)  rJse  the  nor  value  of  £  and  tie  first  sample  set  to  com¬ 
mute  a  new  value  of  00 n  . 

5)  Repeat  steps  3)  and  4)  in  sequence  until  the  successive  values 
of  damping  ratio  and  natural  frequency  either  converge  or 
diverge. 

Although  it  might  seem  that  only  one  sample  set  of  9  and  its  first  two 
derivatives  is  required  to  define  the  parameters  of  a  second  order  system, 
it  is  seen  that  the  method  described  above  requires  two  sets  for  convergence. 
If  only  one  set  were  used  in  both  equations,  the  OJn  obtained  in  the  first 
calculation,  when  substituted  into  the  second  equation  would  produce  the 
originally  assumed  ^  ,  and  convergence  would  be  impossible.  This  is 
because  the  relationships  for  0)n  and  ^  are  actually  just  two  forms  of 
the  same  equation. 

An  appropriate  question  to  ask  at  this  time  is,  will  the  computations 
converge,  and  if  so,  under  what  conditions  will  they  converge?  The  answer 
is  given  in  the  following  section. 
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VI  APPLICATION  OF  PROPOSED  ITS RATI VS  DIGITAL  DAMPING  PATTO  COMPUTER  TO 

A  SECOND  ORDER  SYSTEM 

It  was  felt  that  it  would  be  better  to  test  the  proposed  plant  identi¬ 
fication  method  first  with  a  simple  2nd  order  system  before  adding  the  com¬ 
plexity  of  the  phugoid  and  actuator  modes  of  the  three  degree  of  freedom 
flight  simulation.  When  done  in  this  manner*  the  validity  of  the  basic 
idea  could  either  be  proved  or  disproved  with  certainty  because  the  method 
is  based  on  the  assumption  that  the  airframe  closely  resembles  a  second 
order  system.  Objectives  of  the  program  were  to  determine: 

1)  Will  the  iterative  computing  method  converge  to  correct 
values  of  /  and  cUn  *> 

2)  Conditions  for  convergence, 

3)  Should  both  forms  of  the  equation  be  used,  and  if  so 
when  is  a  particular  form  applicable? 

4)  Effect  of  system  damping  ratio  on  convergence  * 

5)  Effect  of  sample  spacing  on  convergence. 

The  second  order  system  used  in  the  following  programs  was  of  the  form: 

9  +  (DA/VIPiNG  FAC  TOP  )(M)  S  +3G  0  =  0 

where  the  natural  frequency  =  6  radians  per  second  was  approximately  the 
median  of  the  airframe  short  period  CU^  of  this  study,  and  the  damping 
factor  could  be  varied  to  produce  damping  ratios  over  the  range  from  0«1. 
to  i.O  thus  approximating  the  variation  of  the  short-period  mode  damping 
ratio.  Outputs  were  obtained  by  releasing  the  system  from  an  initial  dis¬ 
placement,  0O  at  time  =  0  and  allowing  it  to  oscillate  freely.  Samples  of 

•  m  # 

9,  9j  &  ©  were  periodically  taken  and  stored  for  use  in  the  /  calculations. 
The  dynamic  simulation  was  used  in  four  programs  which  were  performed  to 
test  the  damping  ratio  calculator  and  the  influence  of  certain  variables 
on  its  operation.  The  programs  and  results  are  discussed  in  the  following 
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pages  and  sample  listings  of  each  are  included  in  appendix  VI* 


Program  1.  Calculation  of  U)n  Using  Known  Values  of  £  , 
and  Calculation  of  ^  Using  Known  Value  of  OJn  . 


As  a  check  on  the  accuracy  of  the  dynamic  simulation,  the  samples  of 
output  were  used  in  a  series  of  non-iterative  calculations  in  which  U)n 
was  determined  using  the  known  value  of^?and^  was  determined  using  the 
known  value  of  U)n  . 

It  is  recalled  that  the  relationship  for  natural  frequency, 


G 


Q 


has  two  forms  which  will  be  referred  to  for  identification,  as  the  posi¬ 
tive  and  negative  forms.  Both  forms  were  used  in  the  calculations  in  order 
to  determine  under  what  conditions  each  was  applicable.  A  third  purpose 
of  the  program  was  to  obtain  a  basis  for  comparing  the  accuracy  of  £  and 
U>n  computed  by  the  iterative  method.  Nine  runs  were  made  in  which  damp¬ 
ing  ratio  was  varied  from  0.1  to  0.9  in  0.1  steps,  and  partial  listings 
of  the  runs  made  at  ^  =0.2  and  0.7  are  presented  in  table  IV, 

It  is  seen  from  table  IV  that  the  calculated  values  of  ^  and  CJn 

Q 

have  approximately  0.5%  accuracy  except  in  regions  where  the  ratios  -2-  and 

0 

— g—  become  large.  The  errors  in  ^  and  CJn  obtained  with  this  series  of 
runs  are  due  to  errors  in  the  generated  output  values;  8,  9,  &  0;  and  thus 
are  indicative  of  the  solution  accuracy  of  the  simulated  second  order  system, 
and  they  also  represent  the  highest  accuracy  which  can  be  expected  from  the 
iterative  method.  It  is  also  observed  that  the  calculated  values  of  4  have 
very  little  variation  xchereas  the  u/n  error  increases  in  the  regions  where 
9  is  small,  causing  the  ratios  of  —2—  and  -2—  to  be  large,  and  it  is  in 

C7  C7 

these  regions,  before  9  changes  sign,  that  the  negative  form  of  the  cOn 
equation  is  applicable.  Increasing  damping  ratio  had  the  effect  of  widen¬ 
ing  the  regions  in  which  CJn  errors  appear  and  in  which  the  negative  COn 
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Table  IV  Calculation  of  Natural  Frequency  Given  Correct 
Damping  Ratio,  and  Calculat-ion  of  Damping  Ratio  Given  Correct 
Natural  Frequency  for  a  Second  Order  System  With  CJn-  6. 0  Rad/ 
Sec  and  5=  0.2  and  0.7 

Sample  Spacing  0.0125  Cycle. 


PLANT  IDENTIFICATION  ROUTINE  -  SECOND  ORDER 

CONDI T I 0N=  2  DAMPING  RATI0=  .2000 


ACCELERATION 

VE  L0C I TY 

DISPLACEMENT 

OMEGAN1 

0MEGAN2 

DAMPRAT 

-34.909819 

-.  434  735 

.997322 

6.004203 

-5.829842 

.  1  9.04  94 

-33.595382 

-. 863053 

.989194 

6.004834 

-5.65584 1 

.  1946  19 

-32.  133230 

-1.274004 

.975818 

6.005478 

-5.483248 

.  195985 

-30.53581 1 

-  1. 665819 

.957424 

6.006144 

-5.310185 

.  1966  72 

-28.815956 

-2.036888  .934259 

6.006840 

-5.  134  753 

. 19  7089 

-26.986802 

-2.585762 

.906594 

6.007577 

-4.954950 

.197372 

-25.06 1 70 1 

-2. 71 1 158 

.8747 13 

6.008367 

-4.768574 

.197378 

-23.054149 

-3. 011961 

.838918 

6.009225 

-4.573107 

. 197756 

-20.977702 

-3.  287225 

. 799521 

6.0101 72 

-4.365574 

.  19  7863 

-18.845898 

-3.  536  1  73 

. 756847 

6.01 1234 

-4.  142336 

. 197968 

-16.672187 

-3. 758197 

.71 1229 

6.012447 

-3.8988  10 

. 198057 

-14.469854 

-3. 952858 

.663006 

6.013863 

-3.629053 

. 198134 

-12.251956 

-4.  1  19879 

.612522 

6.015557 

-3.525122 

. 198203 

-10.031248  - 

-4. 259145 

.560125 

6.017644 

-2.9760  75 

. 198264 

-7.820131 

-4.  370  700 

.506159 

6.020513 

-2.566302 

. 198321 

-5.630584 

-4.454739  . 

.450972 

6.023889 

-2.072655 

. 198373 

-3.474 1 19 

-4.511603 

.394904 

6.028996 

-1.459177 

. 198422 

-1.361722 

-4.  541  775 

.330293 

6.036992 

-.666768 

. 198469 

.696184 

-4. 545873 

.281469 

6.051495 

.408725 

.198514 

2.689788 

-4.  524639 

.224752 

6.08634 1 

1.966338 

. 198558 

4.609921 

-4.  478935 

.  168455 

6.272509 

4.362838 

. 198602 

6.448092 

-4.  409  735 

.  1  1 28  7  7 

9 . 794  1  77 

5.852562  - 

. 198645 

..  8.  1 965  1  3_ 

_-4.5  18107_ 

.058305 

23.690226 

5.934103 

. 198688 

1  1.  396610 

-4.072334 

-.046742 

5.972  7  16 

-4  0.82  1967 

.  198  7  78 

12.836410 

-3.  '-20  762 

-.096718 

5.979735 

-22.  194992 

.  198825 

14.  162  73  1 
15.371851 

-3. 751897 

-.  144689 

5.984547 

-16.356640 

.  198875 

-3.567184 

-,  190449 

5.987638 

-13.479787- 

. 198929 

16.459617 

-3.368  1  12 

-. 233809 

5.990128 

r  1  1 .752290 

. 198986 

17.424440 

-3. 156207 

-.274598 

5.992097 

-10.589657 

.199049 

18.264289 

-2. 933022 

-.312667 

5.995  7D8 

-9.745970- 

.199120 

18.978178 

-2. 700125 

-. 347884 

5.995064 

-9.099695 

. 199200 

19.  565846 

-2. 459094 

-.380136 

5.996233  ' 

-8.583825 

. 199292 

20.027742 

-2. 21 1504 

-.409534 

5.997261 

-8. 158339 

. 199402 

20.364998 

-1. 958920 

-.435403 

5.998182 

-7.79782 1 

.199535 

20.579404 

-1.702891 

-.458292 

5.999021 

-7.485314 

.  199704 

20.673579 

-1.444938 

-.477967 

5.999795 

-7.209030 

.199925 

20.649939 

-1. 186546 

-.494414 

6.000520 

-6.960482 

.200234 

2 0.512663 

-. 929164 

-. 507635 

6.001207 

-6.733358 

-200700 

20.265655 

-.  6  74  188 

-.517653 

6. ,001866 

-6.522824 

.201496 

19.913508 

-. 422961 

-. 524506 

6.002506 

-6.325066 

.203192 

19.461264 

-.  1  76  767 

-.528248 

6.003154- 

-6.  136986 

.2094 74 

18.914369 

18.278638 

.063176 

-.52895 1 

6.003  758;- 

-5.955983 

.  168654 

. 295721 

-.526699 

6.004383 

-5.779798 

.  192338 

17.660206 

. 519797 

-.521593 

6.005016 

-5.606394 

.195131 

16.765485 

. 734408 

-.513744 

6.005666 

-5.433857 

.  196223 

15.901  123 

. 938644 

-.503276 

6.006340 

-5.260313 

.  1968  10 

i 
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ACCELERATION 

VELOCITY 

DISPLACEMENT 

0KEGAN1 

_0MEGAN2 

DAM PR AT 

14.973961 

1.  131675 

-.490525 

6.007046 

-5.083842 

.19  7  180 

1  3. 990964 

1.312760 

-.475034 

6.007797 

-4.902394 

.  19  7437 

1  2.959284 

1.481246 

-.457558 

6.008604 

-4.  7  1 369  1 

.  197627 

1  1.8860  12 

1. 636569 

-.438058 

6.009486 

-4.515100 

.  1  97  7  75 

10.778344 

1.7  78253 

-.416701 

6.010463 

-4.303480 

.  19  7895 

9.643431 

1.  9059  14 

-. 393660 

6.01  1564 

-4.074955 

.197795 

8.488369 

2.  01 9255 

-.369  1  13 

6.012829 

-3.824603 

.198080 

7.320155 

2.  1 18068 

-.343239 

6.014315 

-3.545987 

.1981 55 

6. 145654 

2.  202232 

-.316222 

6.016107 

-3.230430 

.  19822  1 

4.971565 

2.  27  1  711 

-.288245 

6.018338 

-2.865862 

.193281 

3.804387 

2. 326550 

-.259490 

6.021224 

-2..434886 

.198336 

2.65059 1 

2. 366876 

-.230141 

6.025 1 54 

-1.911379 

.  198387 

1.515589 

2. 392890 

-.200378 

6.030897 

-1.254149 

. 193436 

.405711 

2. 404869 

-.  170578 

6.040203 

-.594233 

. 198482 

-.673819 

2.403160 

-.  140314 

6.0581 35 

. 792692 

. 198527 

-1.717905 

2.  3881  72 

-.  1  10354 

6.  107544 

2.548848 

.193571 

-2.72 1 799 

2. 360380 

-.030663 

6.565521 

5.  1394  24 

.  1986  14 

-3.681 1 13 

2. 320313 

-.051396 

12.  176150 

5.882218 

.  19865  7 

-4.591836 

2. 268554 

-.022704 

34.023796 

5.944415 

.198701 

-5.450344 

2. 205735 

.005272 

5.964874 

-173.323327 

.198745 

-6.253409 

2.  132527 

.032596 

5.975077 

-32.305450 

.198791 

-6.998207 

2.  049643  j. ,rrit>  .058545 

5.981237 

-  19.985  186 

. 198839 

-7.682317 

1.957825 

.083600 

5.985395 

-15.352948 

.  198690 

-8.303729 

1. 857846 

.  107456 

5.988418 

-12.904145 

. 198945 

-8.860837 

1.  750500 

.130016 

5.990736 

-11.376240 

.  199004 

-9.352439 

1. 636598 

.151191 

5.992589 

-10.322459 

.  199069 

-9. 777735 

1. 516966 

. 170907 

5.994  1  18 

-9.5445 12 

.  199  142 

-10.  136312 

1.392433 

. 189095 

5.995415 

-8.940881 

. 199225 

-10.428141 

1 . 263836 

.205  70  1 

5.996540 

-8.454 1 62 

. 199322 

-10.653365 

1.  132007 

.220678 

5.997534 

-8.049409 

. 199437 

-10.813283 

. 997771 

.233991 

5.998430 

-7. 704088 

. 1995  79 

-10.908339 

.861944  4 

.245615 

5.999248 

-7.402979 

. 19976  1 

-10.940106 

.725327  1 

.255536 

6.000007 

-7.  135387 

.200003 

-10.910265 

>  .588699 

. 263748 

6.000720 

-6.893539 • 

.200347 

-10.820793 

. 452819 

.270257 

6.001598 

-6.67160.3 

.200881 

-10.673937 

.318418 

.275075 

6.002051 

-6.465079 

.201841 

-10.472195 

.  186  199 

.278226 

6.002687 

-6.270382 

.204 1 06 

-10.218301 

.056831 

.279742 

6.003513 

-6.084575 

.216423 

-9.915192 

-. 069053 

.-279662 

6.003936 

-5.905169' 

.  184  186 

-9.565996 

-. 1 90857 

.278033 

6.004563 

%-5. 729981 

.193503 

-9. 174003 

-.308025 

.274909 

6 .005200 

-5.557017  . 

. 195530 

-8. 742644 

-.420043 

.270553 

6.005856 

-5.384384- 

.  196424 

-8.275467 

-.526441 

.264432 

6.006538 

.  -5.210202 

.  196932 

-7.776 1 16 

-.  626795 

.257218 

6.007256 

-5.032525 

.197262 

-7.248304 

-.720726 

.248789 

6.00802 1 

-4.849246 

. 197497 

-6.695796 

-. 807900 

.239228 

6.008848 

-4.658001 

.197673 

-6.  122380 

-.  888033 

.228621 

6.009755 

-4.456031 

. 19 78  12 

-5.531850 

-. 960888 

.217057 

6.010763 

-4.240008 

. 197925 

-4.927986 

-  1. 026274 

.204630 

6.011907 

-4.005796 

.  198020 

-4. 314526 

-1. 084048 

.  191432 

6.013228 

-3.748096 

.1981  02 

-3.695  155 

-1.  1341  12 

.  177561 

6.014792 

-3.4599 1 7 

.  198174 

-3.073483 

-1.  176417 

.163112 

6.016693 

-3. 131756 

.  198239 

-2.453023 

-  1. 2  10955 

.  148182 

6.019084. 

-2. 750264 

. 198297 

-1.837183 

-1.237762 

. 1 32870 

6.0222  18> 

-2.295984 

. 19835 1 

-1.229240 

-1. 256917 

. 117270 

6.026563 

-1.739318 

.198402 

-.632336 

-1. 268539 

. 101478 

6.033073 

-  1.032846 

.  198449 

i 
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CONDITION^  7  DAMPING  RATIO=  .7000 


ACCELERATION 
-32.51 8830 
-29.005397 
-25.687407 
-22.5631 23 
-19.629943 
-16.884500 
-14.322759 
-1 1.940102 
-9. 731413 
-7.691 155 
-5.813440 
-4.092098 
-2.520738 

-  1.092801 

.•198384 
1.359554 
2.397471 
3.318882 
4. 130484 
4.838895 
5.450624 
5.972048 
6. 409392 
6. 763709 
7.055865 
7.276529 
7.436159 
7.539995 
7.593052 
7.600 1 17 
7.565747 
7.494262 
7. 389753 

-  7  Oc;a  >  70 

*  *  j  \j  i  r  " 

7. 096868 
6.915523 
6.715228. 
6.498949 
6. 269440 
6.029253 
5.  78074  1 
5.526069 
5.267215 
5.005985 
4.744016 
4.432787 
4.223622 
3.967704 
3. 716082 
3.469674 
3.229280 
2.995589 
2.769182 
2.550546 
2.340076 
2.  138085 
1 . 944809 


-2. 

-2. 

-2. 

•2, 


__  „  i 


VELOCI  TY 
-.  420371 
-. 804694 

-  1.  1  46321 

-  1. 447687 

-  1. 71  1  196 

-  1. 939218 
-2.  1  34075 
-2.  298034 
-2.433503 
-2. 542022 
-2. 626260 
-2.  6880  10 
-2. 729188 
-2.  75  162  7 
-2.  7570  78 
-2. 747209 
-2. 723603 
-2.  687758 
-2.  64  1  089 

584926 
520520 
449037 
571569 
-2.289  128 
-2.  202652 
-2.  1  13008 
-2.  020993 
- 1. 927337 
•1. 832705 
-1.  737703 
- 1.  642875 
•  1. 548  7  14 
■1. 455656 

uiiaoi 

-  .i  U  ‘t  U  7  ) 

-  1. 274361 

-  1.  186  762 
-1.  101552 

-1.018948 
-.  9  39  133 
-. 862256 
-. 788436 
-. 7 17763 
-. 650302 
-. 586093 
-. 525155 
-. 467489 
-. 413077 
-.361885 
-. 313867 
-. 268962 
-. 227 1 00 
-. 1 88202 
-. 152180 
-. 1 18940 
-. 088382 
-. 060402 
-. 034893 


D  I 


time 


SPLACtM 

.997379 

.989677 

.977440 

.961187 

. 941405 

.918555 

.893063 

.865532 

.835752 

.8046 10 

.772284 

. 73904  7 

. 705169 

.670895 

. 636449 

.602032 

.567826 

.533993 

. 500677 

.468006 

.436089 

.405022 

.374888 

.345754 

.3  176  76 

.29070 1 

.26486 1 

.240183 

.216662 

.  194366 

.  173258 

.  153292 

.  134316 

1  UMOli 

•  I  »  KJ  /  T 

.  1004  06 

.085026 

.070727 

. 057477 

.045242 

.033986 

.023672 

.014262 

.005715 

. 002009 

.00895 1 

.015152 

.020652 

. 025492 

.029712 

.033352 

. 036449 

.039042 

.04 1  166 

.042858 

.044 151 

.045078 

.04567  1 


ENT 


6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 

6. 

6. 

7. 

8. 

9. 

i  n 

I  W  • 

1  1  . 

13. 
15. 
18. 
23. 
29. 
40. 
64  . 
153. 
6. 
6. 
6. 
6. 
6. 
6. 
6. 
6. 
6. 
6. 
6. 
6. 
6. 
6. 


0MEGAN1 
.012666 
.0  1  ?o80 
.012694 
.012708 
.012722 
.0 12737 
.012753 
.012769 
.012786 
.012804 
.012823 
.012845 
.012865 
.012890 
.012917 
.012947 
.01298  1 
.013020 
.015065 
.013120 
.013186 
.013269 
.01 3377 
.013524 
.013737 
.014078 
.014711 
.016515 
.026037 
.5 1 054  1 
.266745 
. 155335 
138630 


?ot;  c;  ao 

^  .J  K/  \  i 

75  70  73 
528653 
792553 
807121 
049256 
507034 
616290 
445604 
293882  • 
012401 
012427 
0.124  50 
0 124  72 
0 12492 
0125  10 
012527 
012544 
012560 
0  125  75 
012590' 
012604 
012618 
012632 


OMEGAN2 
5.422600 
4.874358 
4.370803 
3.904104 
3.467937 
3.057109 
66  7295 
294832 
936570 

t 


6 
6 
6 
6 

'  "  6 

6 

6 

6 

-4  14 

-88 

-4  9 
-54 
-25 
-20 
-17 
-14 
-12 
-1  1 
-9 
-8. 
-7. 
-7 


389752 
.28  1924 
-.92086  1 
-.59450 1 
-.270896 
.051839 
.375869 
. 702179 
1.033626 
1 .3  7  1978 
1  .  7  19973 
2.070580 
2.452076 
2.84  3  1  34 
3.255442 
3.693355 
4.  162070 
4.667840 
5.217939 
5.813242 
6.006  137 
6.0099  10 
6.01 0925 
6.0  1  139  7 

a  m  u7o  • 

I  IV#  I 

6.0  1  185-2 


0  1  1979 
.012075 
.012150 
.012211 
.0  1  226  1 
.012304., 
.0  1 234  1 
.012573 
.409029 
.  150044 
.208001 
.015191 
.886909 
.80  1486 
.302693 
.7554  17 
.761324 
.  187999 
.897905 
815159 
888553 
.082255 


DAMPRAT 

.67  1  594 

.685668 

.690646 

.693036 

.694  'i  78 

.695446 

.696  145 

.696676 

.697096 

.6974 37 

.697721 

.697963 

.698173 

.698558 

.698522 

.698670 

.698805 

.698929 

.699045 

.699153 

.699254 

.699  55  1 

.699444 

.699533 

.6996  19 

.699703 

.699786 

.699867 

.699948 

. 700029 

.700  1  10 

.700192 

.700275 

/flrtU  1 

•  I  *  | 

.  700448 
.  700538 
.700633 
. 70073 1 
.700836 
. 700947 
.701066 
.701 194 
.  70  1355 
.70 1489 
.701662 
.701856 
.  7020  79 
.702538 
.702645 
.705017 
.705479 
.  7040  74 
.  704  8  72 
. 706010 
.  70  7  7  73 
.710899 
.718024 


i 
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\ 


CCEU-RATIO.N 

VELOCITY  DISPLACEMENT 

’  OMEGAN1 

0MEGAN2 

UAMPRAT 

1 . 7604  12 

-.011  745 

-. 045960 

6.0  12645 

-6.3704  15 

. 75 10  1  1 

1.584997 

. 009154 

-.045974 

6.012659 

-5.73589  1 

.637835 

1.418607 

. 027917 

-. 045740 

6.012673 

-5. 158183 

.  6  8  0  o  9  3 

1.261231 

. 044657 

-.045285 

6.012687 

-4.632083 

.  6  8  8  o  0  2 

1 . 1 12814 

. 059486 

-.044632 

6.012701 

-4.  146765 

.691943 

. 973255 

.072514 

-. 043805 

6.012715 

-5.69516 1 

.693794 

.  .842418 

. 085853 

-. 042826 

6.012729 

-3.271520 

.694977 

. 720 133 

. 0936 1 1 

-.041715 

6. .012745 

-2.871087 

.69580? 

.606202 

. 10 1892 

-.040492 

6.012760 

-2.489869 

.6964  12 

.500400 

.  108799 

-.039  1  74 

6.012 777 

-2. 124461 

.696885 

.402483 

. 1 14434 

-.037777 

6.012794 

-1.771913 

.697265 

.312188 

. 1 18895 

-.036318 

6.012812 

-1.429621 

.697577 

.229238 

. 122270 

-.034809 

6.012832 

-1.095246 

.697840 

. 153543 

.  1 24654 

-. 033265 

6.012853 

-. 766647 

.698066 

. 084206 

.  126  132  -.031697 

.  1 26786  time  -.030  1  15 

6.0  128  77 

-.44 1823 

.698263 

.021523 

6.012902 

-. 1 18860 

.698437 

-.035015 

. 1 26695 

-.028530 

6.012950 

.204  1  1  1 

.698594 

-.085719 

. 125935 

-. 02695 1 

6.01296? 

.528959 

.698735 

-. 130900 

. 124575 

-.025384 

6.012998 

.857599 

.698865 

-.  170867 

. 122684 

-.023838 

6.013040 

1.  1920  52 

.698985 

-.205927 

. 120524 

-.022319 

6.013090 

1.534400 

.699096 

-. 236582 

. 1 17555 

-. 020832 

6.013149 

1.88704 1 

.69920 1 

-.262529 

. 1 14432 

-.019382 

6.013222 

2.252559 

.69930 1 

-.284657 

. 1 1 1009 

-.017972 

6.01 35  16 

2.633904 

.699395 

-.303050 

. 107332  , 

-.016608 

6.013440 

3.034477 

.699486 

-.517981 

. 103447 

-.015290 

6.013614 

3.458259 

.6995 74 

-.529718 

. 099395 

-.014022 

6.013876 

3.9099 75 

.699659 

-.538515 

. 095216 

-.012806 

6.014520 

4.395298 

.699742 

-. 344620 

. 090944 

-.01 1642 

6.015248 

4.921011 

.699824 

-.348268 

. 086611 

-.01053? 

6.018406 

5.494219 

.699906 

-. 349685 

. 082246 

-.009477 

6. 158836 

5.9911 20 

.699786 

-.349086 

. 077877 

-.008476 

6.854105 

6.008693 

.700067 

-.346677 

.073527 

-.007530 

7.659854 

6.010503  ' 

.700149 

-.342650 

. 069217 

-. 006638 

8.587596 

6.01 1 184 

.700231 

-. 537189 

. 064966 

-.005799 

9.671900 

6.011543 

.700515 

-.330467 

. 060792 

-.005013 

10.9646  13 

6.01 1765 

.700401 

-.322.644 

. 056709 

-.004279 

12.541722 

6.011917  . 

.700490 

-.513874 

. 052730 

-.003895 

14.52 1378 

6.012028 

.  700582 

-. 504298 

. 048866 

-.002960 

17.097171 

6.0  12113 

. 700679 

-.294047 

. 045126 

-.002373 

20.609855 

6.012180 

.700780 

-.283244 

.041517 

-.001832 

25.7201  14 

6.012236  - 

. 700887 

-.272002 

. 038046 

-.001335 

33.899489 

-  6.012282 

.701002 

-.260425 

. 034718 

-.000880 

49.225428 

6.01232? 

.701  125 

-.248609 

. 031537 

—  . 0004  66 

88.7339  17 

6.012356 

.701259 

-.236640 

. 028304 

-.00009  1 

432.99657-9 

6.012387 

.  7014  06 

-.224599 

. 025621 

. 000247 

6.0-T24  14 

-151 . 100646 

.701568 

-.212556 

. 022889 

.000550 

6.012458 

-64.247932 

.70 1 750 

-.200577 

. 020507 

.000820 

6.012461 

-40.679637 

.701957 

-.  1887  19 

. 0 1 7874 

.001059 

6.012481 

-29.65  1905 

. 702196 

-.  177055 

. 015588 

. 001268 

6.012500 

-23.2297  14 

.702476 

-.  165568 

. 0 1 3447 

. 001449 

6.012518 

-19.006207 

.7028  1  1 

-.  154359 

.01  1448 

. 00 1604 

6.012535 

. -16.002548 

.703221 

-.  143443 

. 009587 

.001736 

6.012552 

9  -13.745592 

. 703739 

-.  132848 

. 007860 

.001845 

6.012567 

-1  1.978633 

. 704418 

-. 122599 

. 006264 

.001933 

6.012582 

-10.550309 

.705354 

-.112717 

. 004794 

.002002 

6.012596 

-9.365561 

.706735 

I 
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equation  is  applicable*  It  would  be  expected  from  the  results  of  this 
series  of  runs  that  the  iterative  method  of  determining  £  and  0Jn  might 

experience  difficulty  in  converging  when  G  is  small  and  that  the  region  of 

non  convergence  will  become  wider  as  system  damping  ratio  is  increasedo 

Program  2.  Iterative  Calculation  of  ^  and  tOn  Using 

Both  Forms  of  the  tdn  Equation. 

The  purpose  of  this  run  was  to  test  the  iterative  damping  ratio  com- 
puter  described  in  section  V.  The  previously  described  second  order  system, 
with  CJn  =•  C  &  =  ,  was  used  to  supply  output  samples  of  9,  8,  and  0 

which  were  stored,  two  sets  at  a  time,  and  used  as  inputs  to  the  comput¬ 
er.  Sample  spacing  in  this  run  was  0.05  cycle  and  each  calculation 

was  iterated  20  times  using  each  of  the  two  forms  of  the  0Jn  equation.  At 

the  termination  of  each  iterating  process,  the  results  of  each  of  the  20 
iterations  was  printed  out  for  study,  the  oldest  set  of  samples  of  9,  0*  & 

9,  discarded,  a  new  set  generated  and  stored,  and  the  calculation  cycle 

repeated  using  the  two  new  sample  sets.  A  sample  series  of  the  printouts 
from  the  £-(A>n  computer  is  listed  in  Table  V  where  the  results  of  one  con¬ 
vergent  and  one  divergent  solution  are  shown* 

The  primary  result  of  this  program  was  the  demonstration  that  the  itera¬ 
tive  method  would  converge,  and  it  did  so  for  51%  of  the  sample  sets  - 
0.2).  Accuracy  obtained  was  of  the  same  order  as  that  obtained  in  Program  #1, 
approximately  0.05%,  and  in  regions  of  small  values  of  6  the  accuracy  of  the 
iterative  method  was  much  greater.  As  an  example,  from  the  iterative 
method  was  6.00359  as  compared  to  =  6.086341  obtained  from  program  #1 
when  the  identical  sample  set  was  used  with  each  method  in  a  region  of 
small  values  of  9.  It  was  found  that  a  necessary  but  not  sufficient  condi¬ 
tion  for  convergence  of  the  iterative  method  was  that  0  and  9  be  of  opposite 


As  is  shown  in  Table  V,  another  interesting  property  of  the  iterative 
method  is  that  when  convergence  occurs,  it  occurs  simultaneously  with  both 
forms  of  the  tOn  equation,  and  the  iterative  series  using  the  negative  form 
of  the  ca)0  equation  converges  to  the  negative  values  of  $  and  arrived 
at  using  the  positive  equation.  Also,  it  was  found  that  the  use  of 
the  negative  ^ n  equation  would  not  be  required  in  future  programs  because 
the  iterations  performed  when  using  the  positive  form  converged  in  areas 

where  the  results  of  Program  #1  indicated  that  the  negative  ^  equation  was 

o 

applicable.  It  had  thus  far  been  established  that*  at  £  «  0.2,  the  itera¬ 
tive  computations  would  converge  for  a  length  of  time  corresponding  to  approxi¬ 
mately  one  quarter  of  the  system  period  at  a  rate  of  two  times  per  system 
cycle.  It  was  felt  that  this  would  suffice  for  the  proposed  adaptive 
system,  however,  it  was  considered  desirable  to  determine  the  effect  of 
different  values  of  damping  ratio  on  convergence. 

Program  3.  Iterative  Calculation  of  /  &  4  at  9 
Values  of  System  Damping  Ratio. 

Nine  runs  were  made  in  this  program,  with  system  damping  ratio  varying 
from  0.1  to  0.9  in  0.1  steps,  and  a  sampling  interval  =  0.05  period.  Only 
the  positive  form  of  the  cJn  equation  was  used  in  the  iterative  process  and 
the  solution  was  tested  for  divergence  after  each  iteration  to  prevent  un¬ 
necessary  calculations.  Table  VI  is  a  partial  listing  of  the  runs  made  with 
/  =  0.2  and  0.7.  It  is  seen  that  the  region  of  divergence  is  appreciably 

wider  with  £  =0.7  than  with  0.2.  The  percentage  of  convergent  iterations 

to  total  iterations  is  tabulated  below  as  a  function  of  damping  ratio  for 
each  of  the  nine  conditions; 
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Ratio .^Uaing^oth^orm^o^the  iJ  Eouatf1  Fre?uen^  and  Damping 
System  «ith  <W>-  6.0™Id/Sec  3  SeC°n<1  0r<3er 

Sample  Spacing  =  0.05  Cycle 


TIME 

0.200003 


-  "f~  Fy-m 


\ 


-Fo 


or /to 


ACCELERATION  VELOCITY  DISPLACEMENT 
-5.63  V  577  -4.454730  0.45097: 


0MEGAN 
7.629993 
-  -  6.818375 
6.413184 
,  6.210091 
6. 1 07890 
6. 056324 
6.0 3020 7 
6.017089 
6.010425 
6.307049 
6.005342 
6.0C4477 
6.D04C40 
6.303819 
6.0;  3707 
6. 3C 3650 
6.303621 
6.0C3607 
6.305599 
6.005595 
0MEGAN 
-2.860373 
-4.239639 
-5.  ''7871  ’ 
-50 528898 
-5. 761842 
-5.880874* 
-5.941390 
-5.  972*‘  06 
-5.987659 
-5.995516 
-5. 999504 
-6.3 J1 522 
-6.002544 
- 6 o  0 .13362 
-6o  3C3323 
-6.003456 
-6.0C3523 
-6.005557 
-6.033574 
-6.003583 


0AMPR 

0. 303380 
0.252439 
0.226C73 
0.212571 
0. 215695 
0.212214 
0. 200434 
0.199537 
0. 199C 83 
0. 198854 
C. 198737 
1.198678 
0. 198648 
0. 198633 
0.198626 
0.198622 
0.198620 
I. 198619 
0.198618 
0.198616 
DAMPR 

0.0761  58 
-0.065534 
-0.132633 
-0. 165552 
-0. 181965 
-0. 19021 y 
r  3 o  194367 
-0. 196468 
-  0  o- 1 9  7  5  3  L 
-0.  1983  67 
r0. 198339 
-0. 198477 
-0.198547 
-0.198582 
-0. 198600 
-0.  198609 
-.'U*  1986  13 
-6*  1 986  16 
-0.198617 
-0. 198617 


9 


/  I 
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Table  V  -  Continued 


f 

TIME  ACCELERATION  VELOCITY  DISPLACEMENT 

0.399999 _ 18.978094 _ -2.700115 _ -0.347882__ 


OMEGAN 

DAMPR 

4.205901 

0.564624 

3.300829 

0.852037 

2.354981 

1.340582 

1.558809 

2.154067 

0.  9881 32 

3.492864 

0.613964 

5.684419 

0.378345 

9.264260 

0.232391 

15.107418 

0.142551 

24.643762 

0.087360 

40.222131 

0. 053459 

65.734811 

0.032669 

107.570.982 

0.019852 

177.027143 

0.011444 

307.084564 

...  _  ..  0.005859 

599.775604 

0.002441 

1439.462204 

-0.002930 

-1199.551773 

44936.139648 

-2894.778473 

108440.639648 

-6985.727661 

261690.726562 

-16858.072266 

OMEGAN 

DAMPR 

-19.190211 

1.053100 

-41.399521 

2.582062 

-97.553423 

6.248351 

-234.412361 

15.085813 

-565.269005 

36.408317 

-1363.942886 

87.862394 

-  -3291.417694 

212.031506  , 

-7942.874573 

511.678204  r 

-19167.864990 

1234.790436 

-46256.204590 

2979.816833 

-111626.239258  ' 

7190.943481 

-269378.281250- 

17353.303223 

-650068.140625 

41877.279297 

-1568755. 15625C 

1C1058. 941406 

-3785745.781250 

243877. 103516 

-9135824. 375000 

588528.257812 

-22046 72C. 500000  • 

1420246.000000 

I487360«  9063S4? 

-128391576.000000 

8270963.6.25000 

...  _  _  -774591 64. 000O0G _ 

19959625,000000 
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Table  VT 


er.eot  Oi  lo'ii  e  Tfe rations  as  a  Function  of  System 

D&r.pjng  ^Mo 


Damping  Ratio 

0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

%  Convergence 

521 

51% 

51% 

49% 

48% 

41% 

34% 

27% 

19% 

The  percentage  of  convergent  iterations  drops  off  sharply  for  values  of 
exceeding  0.5  but  it  is  not  desired  to  increase  the  airframe  damping  ratio 
above  0.7  and  thus  the  minimum  obtainable  ratio  of  convergent  to  total  itera¬ 
tions  should  never  be  less  than  approximately  1/3. 

Program  4.  Influence  of  Sample  Spacing  on  the  Percentage 
of  Convergent  Iterations. 

The  final  program  was  run  to  determine  what  influence*  if  any*  the  width 
of  sample  spacing  has  on  the  ability  of  the  iterative  calculations  to  converge. 
All  runs  were  made  with  £  “0.2,  and  six  values  of  sample  spacing  ranging 
from  0.01  cycle  to  0.5  cycle  were  investigated.  Results  of  the  program  are 
tabulated  below. 

Table  VIII  -  Percent  of  Convergent  Iterative  Calculations  as  a  Function  of 
Sample  Spacing  /  =0.2 


Spacing 

0.01  cycle 

0.02 

0.05 

0.1 

0.2 

0.5 

%  Convergence 

56% 

53% 

51% 

50% 

48% 

46% 

It  is  seen  that  narrow  sample  spacing  favors  convergence 9  however  it  can  be 
varied  over  wide  limits  with  little  effect  on  convergence.  It  was  arbitrar¬ 
ily  decided  to  set  sample  spacing  at  0.025  seconds  for  the  three  degree  of 
freedom  airframe  simulation;  the  primary  consideration  being  conservation  of 
computer  time.  This  sample  spacing  corresponded  to  approximately  50  samples 
per  short  period  cycle,  depending  on  flight  condition. 
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Table  VII  Iterative  Calculation  of  Natural  Frequency  and  Damping 
Ratio  of  a  Second  Order  System  With  ^=6.0  Rad/Sec  and  /=  0.2  &  0.7 
Sample  Spacing  =  0.05  Cycle 


Condition  2 

Damping  Ratio  « 

0.2 

TINE 

ACCELERATION 

VELOCITY 

DISPLACEMENT 

0.2C0000 

-5.630577 

-4.454730 

0.450971 

ONEGA-N 

DAMPRAT 

7.629993 

0.303380 

6.818375 

0.252439 

6.413184 

0.226073 

6.210091 

0.212571 

6.107890 

0.205695 

6.056324 

0.202204 

6.030267 

0.200434 

6.017089 

0.199537. 

6.010423 

0.199083 

6.007049 

0.198854 

6.005342 

0.198737 

6.004477 

0.198678 

6.004040 

0.198648 

6.003819 

0.198633 

6.003707 

0.198626 

6.003650 

0.198622 

6.003621 

0.198620 

6.003607 

0.198619 

6.003599 

0.198618 

6.003595  __  C. 198618 


TIME  ' 

ACCELERATION 

VELOCITY  DISPLACEMENT 

0.250000 

2.689775 

-4.524627  0.224752  .  • 

0MEGA-N 

DAMPRAT 

9.251958 

0.261913 

6.966586 

0.215692 

*«  ^ 

6.256753 

0.202902 

6.066626 

0.199669 

6.019025 

0.198874 

*• 

6.007354 

0.198680 

6.004508 

0.198633 _ 

- 

6.003815 

0.198622 

6.003646 

0. 19861 9_ 

6.003605 

0.198618 

6.003595 

0.198618 

6.003592 

0.198618 

6.003592 

0.198618 

.  Jr 

6.003592 

0.198618 

6.003592 

0. 198618_ 

6.0C3592 

0.198618 

6.003592 

r  0.  198618  _ 

6.003592 

*  0.198618 

i 
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Table  VII-  Continued 


Damping  Ratio  =  0.2 


6.003592  0.198618 

6.003592 _  0.198618 


TIME 

ACCELERATION 

VELOCITY 

DISPLACEMENT 

0.299999 

9.848088 

-4.205210 

0.005013 

OMEGA-N 

DAMPRAT 

15.324385 

0.085544  _ 

-7.279872 

-0.165185 

-4.234596 

_  -0.279041 

-1.191589 

-0.983381 

-SOLUTION  DIVERGING 


TIME 

ACCELERATION 

VELOCITY 

DISPLACEMENT 

0.349999 

1 5.  371489 

-3.567171 

-0.190448 

OMEGA-N 

DAMPRAT 

SOLUTION  DIVERGING 

TIME 

ACCELERATION 

VELOCITY 

DISPLACEMENT 

0.399999 

_ L8.-978094 

-2.7QQ 11 5 

-0.347882 

OMEGA-N  DAMPRAT 

_ 4.205901 _ ; _ 0.564624 

3.300829  0.852037 

-SOLUTION  DIVERGING  _ 


TIME 

ACCELERATION  VELOCITY 

DISPLACEMENT 

0.449999 

-  20.579304  -1.702884 

-0.458290 

OMEGA-N 

DAMPRAT 

_ 4.907358 

0.570963 

4.181908 

0.882182 

• 

SOLUTION  -DIVERGING- - 

TIME 

ACCELERATION 

VELOCITY 

DISPLACEMENT 

0.499998 

__20. 265548 _ 

-0.674185 

-0.517650 

OMEGA-N 

DAMPRAT 

"  5.377643 

0.130322 _ 

SOLUTION  DIVERGING 


TIME _ ACC  EL  E  RAT  I  ON VE  L  OC I T  Y _ D I S  P  LACEMENT_ 


0.549998 

18.278530 

0.295719 

-0.526696 

...  QM.EGA-N 

DAMPRAT 

4 

5.757618  -0.240352 

SOLUTION  DIVERGING _ 


i 
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Table  VII-  Continued 


Damping  Ratio 


0.2 


TINE 

ACCELERATION 

VELOCITY  DISPLACEMENT 

0.599997 

14.973863  _ 

1.131666  _ -0.490321 _ 

QMEGA-N 
6.11 9884 

DAMPRAT 

0.244751 

6.030041 

0.209181 

6.009638 

0.201036 

6.004975 

0.199172 

6.003908 

0.198745 

6.003664 

0.198647 

6.003608 

0.198625 

6.003595 

0.198619 

6.003593 

0.198618 

6.003592 

0.198618 

6.003592 

0.198618 

6.003592 

0.198618 

. 

6.003592 

0.198618 

6.003592 

0.198618 

6.003592 

0.198618 

6.003592 

0.198618 

6.003592 

0.198618 

6.003592 

0.198618 

6.003592 

0.198618 

6.003592 

0.198618 

TIME 

ACCELERATION 

VELOCITY 

DISPLACEMENT 

0.649996 

10.778271 

1.778240 

-0.41  6698___ 

ONEGA-N 

DAMPRAT 

6.525987 

0.300233 

6.262415 

0.249806 

6.132749 

0.224382 

6.068288 

G. 211580 

6.036061 

0.205138  _ 

6.019904 

0.20189.7 

..  -  " 

6.011790 

0.200267 .  _ 

6.007714 

0.199447 

- 

6.005664 

0.199035 

6.004634 

C. 198828 

6.004116 

0.198723 

6.003855 

0.198671 

.  * 

6.003724  .. 

G.  1 98645 

6.003658 

C. 198631 

6.003625 

G.  198625  __ 

_ __ _ . 

6.003608 

0.198621 

6.003600 

0.198620 

6.003596 

0.198619 

.  jr  - 

6.003594 

0.198618 

6.003593 

0.198618 

i 
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Table  VIE-  Continued 

Condition  7  Damping  Ratio  =  0.7 


TINE 

„ACCELERAT  ION  VELOCITY 
'  1.  359547  _  -2.747204 

DISPLACEMENT 

0. 200000 

0.602031 

ONEGA-N 

DAKPRAT 

4.221364 

_ 0.  521.158  _ 

4.917109 

0.589098 

5.325055 

0.629943 

5.575466 

0.655293 

5.732646 

0.671299 

5.832533 

C.  68 1505 

5.896477 

0.  688051 

5.937594 

0.692266 

_5_.-9_6.41 07 

0.694986 

5.981233 

C. 696743 

5.992308 

.  C.  69788C  ... 

5.999476 

0.698616 

6.004117 

0.699093 

6.007122 

0.699402 

6.009069 

0.699602  . . 

6.010331 

C. 699731 

6.011  148 

0.699815 

6.011 677 

0.699870 

6.012021 

€.699905 

6.012243 

C. 699928 

TINE 

ACCELERATION  VELOCITY 

DISPLACEMENT 

0.250C00 

4.838883  -2.584920 

0.468005 

ONEGA-N 

DANPRAT 

2.861347 

0.586139 

4.887297 

P. 633940 

5.364674 

0.660114 

5.622871 

C. 675476  . 

5.773553 

C. 684772 

t 

5.864458 

0.690488  _ 

5.920257 

0.694035 

_5. 954843 

C. 696248 

5.976406 

0.697633 

5.989896 

0.69850 1  . 

5.998354 

0.699046’ 

6.003663 

0.699389 

6.007000 

0.699605 

6.0C9C97_ 

0.6  99740  .  .. 

6.010416 

0.699825 

• 

6.011245 

0.699879 

6.011767 

0.699913 

6.012096 

_0. 699934 
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Table  VII-  Continued 


Damping  Ratio 


0.7 


6.012302 

6.012432 


u. 699947 
0.699956 


TINE 

ACCELERATION 

VELOCITY 

DISPLACEMENT 

0.299999  6.768687 

OMEGA-N  CAM  PRAT 

-3.249:22  - 0 . 7 ' 0  4 1 2 

-1.717651  -0.990456 

SOLUTION  CIVERGING 

-'>.289121 

• 

0.345753 

TIME 

0.349999 

OMEGA-N 

-9.915033 

-2.306982 

ACCELERATION 
7.539968 
CAM.PRAT 
-0.  81  5>.'83 
-C. 991656 

VELOC I TY 
-1 .927330 

DISPLACEMENT 

0.240162 

SOLUTION  CIVERGING 

TIME 

0.399999 
OM.EGA-N  ....  . 

ACCELERATION 

7.494233 

CAM.PRAT 

VELOCITY 

-1.548708 

DISPLACEMENT 
j. 153291 

SOLUTION'  DIVERGING 

_  TIME  _ 

ACCELERATION 

VELOCITY 

DISPLACEMENT 

0.449999 

ONFGA-N 

6.915493 

CAMPRAT 

-1  .  186757 

0.085026 

SOLUTION  CIVERGING 

TIME _ 

A  CC  EL  ER AT  I  ON  V  EL .00  I  T  Y 

DISPLACEMENT 

0.499998 

OMEGA-N 

6.029223 

C  A  M  P  R  A  T 

-P.862252 

0.033986  7 

-44.580629 
SOLUTION  CIV5RG 

-0.957002 

INC-.  _ 

- 

-• 

TIME 

.0.549998 

ACCELERATION 
5. Of  5958 

VELOC  IT Y 
-P .586093 

DISPLACEMENT 

•  '  -0,002009 

CMEGA-N 

-64.266227 
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G.  04370.1  _ 
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_ TIME 

ACC.ELEE.At  ION 
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0.599997 
OMEGA-N  •, 
1C. 228298 
20.289338 

3.967681 

CAMPRAT 

0. 17571 1 
-  .444423  ._ 

-''.361  883 

-7.C25492 

SOLUTION  CIVERGING 
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Table  VII-  Continued  Damping  Ratio  »  0.7  * 


_ TINE 

0.649996 

OVEGA-N 

ACCELERATION. 

2.995569 

C.AMPRAT 

0.158457 

-0.318348 

.VELOCITY _ 

.  1 8820C 

DISPLACEMENT 
-0.03904 1 

8.028860 
_  10.427491 

SOLUTICN  DIVE 

RGING 

TINE- 

ACCELERAT ION 

VELOCITY 

DISPLACEMENT 

0.699996 

2.  138C70 

.060402 

-0.045078 

OMEGA-N 

DANPRAT 

7.040965 

-C. 113642 

SOLUTICN  DIVERGING 

TIME 

ACCELERATION  VELOCITY 

DISPLACEMENT 

0.749995 

1.418596  ".G27917 

-0.045740 

_ ONEGA-N 

DANPRAT 

SOLUTICN  DIVERGING 

TINE 

ACCELERATION  VELOCITY 

DISPLACEMENT 

C. 799994 

0.84241  1  ''.083853 

-0.042826 

CNEGA-N 

CAMPRAT 

• 

5.818552 

0.622535 

5.961974 

0.679927 

0 

5.999497 

0.694778 

6.0C9243 

0.698625 

6.011769 

0.699621 

_ 6.012424 

0.699880 

6.012594 

0.699946 

6.012637 

_jQ,  699964 

6.012649 

0. 699966 

T* 

6.012652 

0.699969 

6.012653 

0.69997C 

6.012653 

0.699970 

6.012653 

C. 699970 

_ 6.012653 

0.699970- 

6.012653 

0.69997 T 

6.012653 

_  0.699970 _ ...  ...  .....  .  . 

6.012653 

0.699970 

6 . C 1 2  653 

0.69997  0 

6.012653 

£>.□ 

0.699970 
...  J8.69.99M_ 
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VII  -  APPLICATION  OF  THE  DAMPING  RATIO  COf II  UTER  TO  A  SELF  ADAPTIVE  PITCH 

DAMPING  SYSTEM 

Investigation  of  the  Ability  of  the  Computer  to  Determine 

Short-Period  Mode  Characteristics 

It  has  thus  far  been  established  that  the  damping  ratio  computer  opei- 
ates  satisfactorily  when  supplied  input  samples  generated  by  a  second  order 
system,  and  it  has  also  been  established  that  the  short  period  mode,  a 
second  order  oscillation,  is  the  dominant  mode  of  the  airplane  longitudinal 
dynamics.  The  next  step  in  the  program  is  then  to  investigate  the  opera- 
tion  of  the  computer  when  sut  plied  samples  generated  by  the  three  degree  of 
freedom,  fourth  order  airframe  simulation.  The  simulation  program  is  dis¬ 
cussed  in  detail  in  Appendix  VI. 

Flight  condition  number  seven # M  =  1.8  ?  25,000  ft.,  having  short-period 
mode  damping  ratio  =  0.095  and  period  =  1.2  seconds,  was  chosen  for  the  test. 
At  the  start  of  the  computer  run,  the  airframe  was  released  from  a  non-trim 
pitch  condition  and  alio wed  to  oscillate.  Periodic  samples  of  pitch  angle, 

0,  and  Its  first  two  derivatives  were  taken  and  supplied  to  the  com¬ 

puter  as  described  in  the  preceeding  section.  At  the  beginning  of  each  set 
of  calculation  an  initial  value  of  f  =  0.4  was  introduced  into  the 

equation  and  the  iterations  were  carried  out  15  times.  The  iterations  con¬ 
verged  for  49  percent  of  the  sample  sets,  which  is  in  good  agreement  with 
the  previous  second  order  system  investigation,  however  the  computed  values 
of  damping  ratio  ranged  from  0.1448  to  0,024565  during  a  2.5  second  (flight 
time)  run.  Such  a  wide  variation  of  computed  £  was  inconsistent  with  the 
results  of  the  preceeding  study  in  which  it  was  demonstrated  that  the  longi¬ 
tudinal  aircraft  dynamics  very  closely  approximate  the  dynamics  of  a  second 
order  system,  A  review  of  foregoing  investigations  was  conducted  to  deter¬ 
mine  if  there  were  any  significant  differences  between  the  second  order 
study  and  the  present  application  of  the  computer,,  and  the  answer  was 
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found  in  c xamtn -t 5.  >n  <  *  airrlar.?  short-period  titrh  angle- to-e levator 

deflection  transfer  tuner 1 o 

J9  =  (S+Te) _ 

It  is  seen  that  the  l/S  term  causes  pitch  angle  output  to  exhibit  a  theo¬ 
retically  infinite  steady  state  error  or,  if  the  ait  frame  is  perturbed  from 
its  trim  position  at  a  given  Flight  condition*  the  ensuing  short  period  oscil 
lations  will  not  decay  to  the  original  trim  angle.  Further*  it  can  be  shown 
that  the  airframe  will  return  to  trim  condition  as  a  function  of  the  longi¬ 
tudinal  phugoid  mode.  The  primary  difference  between  the  ideal  second  order 
system  and  the  short  period  dynamics  is  then  the  fact  that  the  former  system 
oscillates  about  0  equal  zero  whereas  the  airframe  oscillates  about  an  aver¬ 
age  pitch  angle  which  is  constantly  varying  as  a  function  of  the  phugoid 
dynamics,  and  thus  the  pitch  angle  samples  introduced  into  the  computer  were 
erroneous.  The  solution  to  this  problem  was  seen  in  the  pitch-rat e-per-stabi 
lizer-def lection  transfer  funct ion, 

Q  _  _ tie  (£  Ts ) _ 

lt  (6*--h2.fU)nS>-h  U)nz) 

which  exhibits  zero  steady  state  error,  and  which  also  shows  that  free  air¬ 
frame  short  period  oscillations  oscillate  about,  and  decay  to,  zero  pitch 
ratet  It  was  realized  that  pitch  rate  is  the  true  second  order  variable  of 

the  airframe  longitudinal  dynamics  and  in  the  study  to  be  outlined  next, 

•  *  •  •  • 

pitch  rate  and  its  first  two  derivatives  9  and  9  were  utilized  as  inputs  to 
the  %-uJn  computer, 

A  second  test  of  the  $-u)n  computer  using  the  three  degree  of  freedom 
flight  simulation  was  performed,  as  stated  above,  and  a  partial  listing  of 
the  results  obtained  when  sampled  values  of  8,  0  &  0,  were  used  as  computer 
inputs  is  presented  in  Table  IX.  As  before,  the  calculations  converged  for 
49%  of  the  sample  sets  but,  as  can  be  seen,  the  calculated  values  of  damp¬ 
ing  ratio  compare  favorably  with  the  known  £  =  0.095.  It  can  also  be  seen 
from  examination  of  Table  IX  that,  in  most  cases,  the  calculations  have  not 

7, 


quite  settled  out  to  sre«..dy  \  <)  c  £  and  CUn  aft  15  i  ter&t  lonr,^  and 

as  a  result,  It  was  decided  to  re* eat  the  iterative  calculation  cycle  20 
times  in  all  future  programs* 

Application  of  the  computer  to  the  Adaptation  Loop 

It  was  felt*  at  this  time,  that  the  computer  was  generating  suf¬ 

ficiently  accurate  values  of  damping  ratio  to  allow  its  use  in  the  adapta¬ 
tion  loop,  however  one  obvious  discrepancy  remained  to  be  corrected  before 
the  loop  could  be  closed*  That  is*  what  should  be  clone  regarding  the 
calculations  which  do  not  converge  and  thereby  generate  highly  inaccurate 
values  of  (  ?  Before  discussing  the  technique  which  \/as  used  to  detect  and 
account  for  divergent  cal  cl lat ions*  possibly  a  better  understanding  may  be 
obtained  from  a  short  review  of  the  proposed  mechanization  of  the  adapta¬ 
tion  loop* 

It  will  be  recalled  that  it  was  proposed  to  compare  the  output  of  the 
damping  ratio  computer  with  the  desired  value  of  short  period  damping  ratio 
*  the  error  between  £  -  desired  and  £  -  computer  being  used  as  the  input 
signal  to  a  rate  servo  which  would  vary  pitch  damper  gain  in  the  proper 
direction  to  null  out  the  error,  and  at  a  rate  proportional  to  the 
magnitude  of  the  error*  It  is  apparent  that  the  erroneous  values  of  £ 
resulting  from  divergent  calculations  could  not  be  allowed  to  reach  the 
computer  output,  and  it  was  found  that  this  could  be  accomplished  by 
simply  including  a  routine  which  examined  the  resultant  values  of  Ldn  and  £ 
at  the  termination  of  the  series  of  iterative  calculations  for  each  sample 
set.  After  examination  of  a  large  number  of  past  output  records  from  both* 
the  second  order,  and  the  airframe  studies,  it  was  determined  that,  in 
every  case,  when  the  calculations  diverged  the  resultant  calculated  values 
of  either  £  or  COn  or  both  were  completely  out  of  the  range  of  variation 
of  the  real  values  of  £  and  U)n  *  For  examples  of  divergent  calculations, 
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Table  IX  Iterative  Calculation  of  Airframe  Short  Period  Damping 

?ti^in«?d«piatUrS1  Fre^uln°y  Using  ©,  0,  &  0,  Generated  by  a  iJng- 

Inputs  1  T  rSe  egree  °f  Freedora  Flight  Simulation,  as  Computer 


Natural  Frequency  =  5.24  Rad/Sec 
Damping  Ratio  »  0.095 


TINE-SEC 

THETA 

THETACOT  F4D  VELOCITY 

.  —1 «  09599_9_  . 

0.  744P06 
CNEGAN 

7o 98.18' 

1  0  7  7  6  2  2  4  80  5 .  9  8213  46,  . 

DA~MPR 

,  #  ~  6  Y  7 ,4  9 

70  226988 
.  6o.aE.5L4 

. '0  266  7  ‘ 

•  .  7  7  r  P  7 

6,294  84 

60  ’  6*4 

*  1898  97 

. .  0  I  6  5376 

j  a  1  4  7  :  6 

5  0  8  4166 

5  0  6  5  3 1 r; 

So  542262 
.5 o-4b  £81 

—  .0  1  J  O  .  C  1 

>  10  7,  7  5 

u  7  1  ^  <  7>< 

- - - -  _ 

5o4  26 

5 0  3.35  62  7 

5  ©  3  2  2  5  3 '» 
3o998  * 

-  - - - — — — -  C  i.  I.-.-  9__ 

0  i  96  69 
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'.1 .  1  7  1  8  c  p 

60  '67493 
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Table  IX 


Continued 
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see  Tables  V  or  VII.  a1 so,  t  was  known  with  certainty,  that  for  the  air¬ 
craft  considered  in  this  study,  damping  ratio  was  restricted  well  within 
the  range  from  t\  to  1,0  and  natural  frequency  within  the  range  from  0* 
to  8.0  radians/second.  Thus  a  divergence  detection  routine  was  programmed 
to  use  a  simple  logical  process  for  determining  whether  or  not  the  calculated 
values  of  both  and  CUn  fell  within  the  allowable  ranges.  If  both  values 
were  acceptable  to  the  detection  routine,  the  calculated  value  of  £  was 
set  into  the  damping  ratio  comparison  process,  and  if  either  £  or  OJn  or 
both  were  unacceptable  to  the  detection  routine,  the  last  acceptable  com¬ 
puted  value  of  %  was  reset  into  the  comparator.  The  values  of  U)n  and  £ 
resulting  from  a  divergent  series  were  always  orders  of  magnitude  different 
from  the  allowable  ranges  of  u)n  and  £  ,  and  thus  the  allowable  ranges 

could  be  altered  considerably  to  suit  the  characteristics  of  other  systems 
with  no  deterioration  in  the  effectiveness  of  the  divergence  detection 
routine.  As  an  example,  if  it  were  expected  that  the  airframe  might  have 
short-period  poles  in  the  right  half  plane  over  some  flight  conditions,  the 
allowable  £  range  could  be  increased  from  £  =  -1.0  to  +1,0. 

A  variable  which  has  not  yet  been  discussed,  but  which  was  quite  im¬ 
portant  in  the  evaluation  of  the  adaptation  loop,  is  the  gain  of  the  rate 
servo  (hereinafter  called  the  adaptation  servo)  used  to  vary  pitch  damper 
gain.  When  ^ -  error  signals  are  supplied  to  the  adaptation  servo  which 
in  turn  causes  variations  in  airframe  damping  ratio,  the  entire  adaptation 
loop  appears  to  resemble  a  second  order  servomechanism  and  should  be  amenable 
to  conventional  systems  analysis  techniques.  However  analysis  was  compli¬ 
cated  by  the  intermittent  action  of  the  damping  ratio  computer  and  thus  an 
"optimum”  value  of  adaptat ion-sei vo-gain,  Kg,  was  determined  experimentally 
by  performing  a  series  of  three  degree  of  freedom  computer  runs  with  vary¬ 
ing  values  of  K  -  Flight  condition  7,  M  =  1,8  »? . 25,000  ft,,  was  again  used 
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Figure  29.  Time  History  of  Airplane  &  Control  Surface  Motion 
and  Pitch  Damper  Gain  During  Adaptation  Process  With  No  Adapt¬ 
ation  Servo  Lock 

Initial  Kfl  =  0,  Adaptation  servo  gain  =  1. 


i 
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for  c.  1 1  runs  and  adaptation-servo  gain  was  varied  fcom  K  «  0.05  to  2.0 

s 

in  six  steps. 

Figure  29  is  a  typical  time  history  of  the  run  made  with  K  *  1.0. 

s 

Airplane  pitch  angle,  pitch  rate*  control  surface  position*  and  pitch  damp¬ 
er  gain  are  shown  as  functions  of  time*  and  again  the  run  was  started  by  re¬ 
leasing  the  airframe  from  a  non-trim  p>itch  position  and  allowing  it  to  oscil 
late.  The  damping  ratio  of  the  basic  airplane  at  this  condition  was  0.096 
and  the  steady  state  damper  gain  for  correct  adaptation*  from  the  root  locus 
plot,  was  u.344.  Initial  damper  gain  was  set  at  zero  at  time  =  zero  and  it 

is  seen  that  K,  =  0.344  was  achieved  in  0*45  seconds*  and  that  the  airframe 
a 

motion  xca s  heavily  damped  by  the  control  action.  However,  it  is  also  seen 

that  damper  gain  did  not  level  out  at  the  desired  value  of  R.  =  0*344  but  in 

d 

stead  continued  increasing  in  a  somewhat  random  manner  until  at  time  ~  10 
seconds  (not  shown)  it  reached  the  required  value  for  instability  of  the 
actuator  mode  and  a  violent  divergent  oscillation  ensued.  In  this  condition 
the  actuator  motions  comprised  the  dominant  longitudinal  mode,  and  it  was 
interesting  to  note  that  the  computer  was  still  performing  its  duty  in 

that  it  was  now  computing  natural  frequency  and  damping  ratio  of  the  actua¬ 
tor  mode.  The  crossover  frequency  for  this  mode  was  predicted  to  be  34 
radians  per  second  from  the  root  locus  plot,  which  was  verified  by  the  fact 
that  computed  damping  ratio  changed  from  positive  to  negative  sign  when 
computed  COn  was  approximately  34,5  rad/ sec.  But,  the  operation  of  the 
adaptation  loop  was  obviously  unacceptable  with  =  1.0.  An  improvement 
in  performance  was  detected  when  lower  values  of  adaptatlon-servo-gain  were 
investigated,  and  with  1C,  ~  0.1  the  damper  gain  did  settle  out  reasonably 

a 

well,  but  approximately  eight  seconds  were  required  to  complete  the  adapta¬ 
tion  process  when  using  the  K  .  Osburn  has  shown  in  Ref.  127  that  the  MIT 
self  adaptive  system  is  capable  of  completing  the  adaptation  process  on  a 
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similar  airfr;  e  luring  l v, o  cycles  or  \  lam  :)Sci1Uqo  which  would 
cone">[Oncl  to  adaptation  tine  -  2,4  seconds  for  flight  condition  7, 
presently  under  discussion*  Thus,  a  review  was  made  of  the  adaptation 
method  to  determine  if  the  process  could  be  speeded  up  without  attend¬ 
ant  drift  in  pitch  damper  gain* 

The  inability  of  the  adaptation  loop  to  set  proper  damper  gain  when 
operating  at  high  adar  tat  ion-servo-gains  was  fount!  to  be  caused  by  four 
interact ing  factors : 

1}  The  reader  will  recall  that  when  the  process  failed  to 

converge,  the  computer  set  the  most  recently  computed  con¬ 
vergent  value  of  fy  into  the  damping  ratio  comparator,  and 
the  adaptation  servo  continued  co  change  damper  gain  at  a 
rate  proportional  to  the  indicated  £  -  error* 

2)  It  was  found  that  during  the  adaptation  process,  the  const¬ 
antly  changing  damper  gain  caused  errors  in  computed 
and  £  due  to  the  fact  that  system  parameters  were  altered 
appreciably  during  one  sampling  interval.  And.  the  comput¬ 
ed  values  of  4  were  caused  to  vary  in  a  random  manner  and 
were  generally  0*C5  to  0*1  below  the  true  value  of  damping 
ratio, 

3}  As  pitch  damper  gain  increased,  the  effect  of  the  time  rate 
of  change  of  system  parameters  on  the  calculations 

was  magnified  due  to  the  increased  system  loop  gain, 

4)  when  system  motions  damped  to  zero,  the  computer  was 

not  supplied  with  sufficient  information  to  perform  calcu- 
lat ions , 

It  can  be  seen  how  these  factors  combined  to  cause  erratic  operation  of  the 
adaptation  mechanism.  As  an  example,  if  a  value  of  4  =  C,6  were  computed 
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and  set  intc  the  £  comparator,  and  no  new  convergent  values  of  £  were 
geneixited  for  two  seconds*  a  steady  £  -  error  signal  =  0,1  would  be  fed 
to  the  adaptation  servo  and  it  would  increase  damper  gain  at  a  constant 
rate  during  the  two  second  period.  This  is  Illustrated  in  Figure  29*  where 
was  Increased  linearly  from  time  »  3.5  seconds  to  5*5  seconds 0  At  5.5 
seconds  a  convergent  value  of  £  ,  exceeding  0*7,  was  obtained,  set  into  the 
£  -  comparator,  and  damper  gain  began  a  downward  drift*  However*  the  drift 
of  Kd  was  generally  upward  because  the  erroneous  damping  ratio  solutions 
were  usually  smaller  than  0.7, 

The  problem  of  drifting  damper  gain  was  solved  by  an  alteration  of  the 
computer  output  routine  in  which  output  logic  was  re-programmed  to 
set  $  =  0,7  into  the  f  comparator  if  no  convergent  solution  were  obtained. 
This  action  effectively  "locked"  the  adaptation  servo  at  its  present  set¬ 
ting,  and  damper  gain  was  not  allowed  to  vary  until  generation  of  a  new 
convergent  solution.  The  above  process  will  be  referred  to  as  the  "adapta¬ 
tion-servo  lock"  in  the  remainder  of  this  report.  An  important  benefit 
accruing  from  the  use  of  the  servo  lock  was  the  fact  that  the  plant  was 
allowed  to  settle  to  a  linear  system  during  each  period  of  nonconvergent 
calculations  and  thus  the  accuracy  of  subsequent  convergent  calculations 
was  improved.  It  is  recalled  that,  when  the  plant  is  in  motion*  the  calcu¬ 
lations  converge  during  one  quarter  period  in  each  half  system  cycle*  and 
from  this,  it  would  be  expected  that  damper  gain  T7ould  be  changed  in  quart¬ 
er  period  steps  during  the  adaptation  process.  This  is  borne  out  in  Figure 
30  I'/hich  is  a  time  history  of  the  adaptation  process  with  the  servo  lock  In 
operation. 

The  run  shown  in  Figure  30  is  identical  in  all  other  respects*  except 
the  use  of  the  servo  lock,  to  the  run  illustrated  in  Figure  29  that  is, 


32 


ni«ht  con:’ it  Lon  7  ,  .uk1  n-'a\  ti«.  ion-servo  3a  in  1,./.  Tt  is  saen  that 
pitch  damper  gain  increased  in  steps,  as  predicted,  and  it  settled  out 
close  to  the  predicted  valae  o'"  K.  =  0,3V*  in  2.3  seconds.  After  three 

vi 

seconds  running  time,  the  airplane  motion  damped  cut  and  no  values  of  £ 
were  available  until  t  =  5  sec.  During  this  interval  damrer  gain  remained 
constant  at  Y  p  =  0.?*  At  t  -  5  sec.,  the  airframe  was  perturbed  by  an  im¬ 
pulse  and  the  computer  begcan  generating  values  causing  damper  gain  to  be 
moved  into  the  Y ,  ~  0.34  to  0,35  region  again.  The  time  history  of  notions 

Ci 

resulting  from  the  second  impulse  also  serves  to  illustrate  the  behavior  of 
the  adapted  airframe  which  displays  a  dariping  ratio  =  0.7. 

A  series  of  runs  was  male  with  the  servo  lock  in  operation  using  eight 
different  values  of  adaptation  servo  gain  ranging  front  Ks~  0.1  to  5.0.  It 
was  found  from  this  series  that  values  of  Tf  exceeding  1.0  caused  an  ex- 

o 

cessive  overshoot  in  damper  gain  during  the  adaptation  process  and  that 
values  of  K  less  than  1.0  resulted  in  a  smooth  adaptation  transition,  but 
time  required  to  fully  adapt  was  approximately  inversely  proportional  to 
servo  gain.  On  the  basis  of  these  results,  it  was  decided  that  Kg  ~  1.0 
was  optimum,  and  a  second  series  of  runs  was  conducted  to  investigate  the 
operation  of  the  adaptation  loop,  with  this  gain  setting,  at  each  of  the 
11  flight  conditions.  Table  X  is  a  sample  listing  of  the  computer  output  for 
the  three  degrees  of  freedom  run  performed  for  flight  condition  7,  from  which 
the  reader  can  see  the  stepwise  manner  in  which  adaptation  is  accomplished. 
Items  of  interest  in  this  program  were  the  time  required  to  adapt  the  air¬ 
frame  from  its  basic  un-nugrr.ented  condition  to  damping  ratio  =  0.7,  and  the 
steady  state  value  of  damper  gain  upon  completion  of  adaptation.  These  values 
are  listed  in  Table  XI  for  each  flight  condition.  Other  Items  listed  are  the 
theoretically  pre tinted  gain  for  adaptation  and  the  basic  airframe  short  per¬ 
iod  damping  ratio.  Fgiure  30  shows  that  the  steady  state  damper  gain  exhibit¬ 
ed  small  amplitude  fluctuations  about  a  mean  value  during  operation^  and  it  Is 


Figure  30.  Time  History  of  Airplane  &  Control  Surface  Motion 
and  Pitch  Damper  Gain  During  Adaptation  Process  With  Adapt¬ 
ation  Servo  Lock 

Initial  K^=  0,  Adaptation  Servo  Gain  -  1. 
Flight  Condition  7  -  M=  1.8  @  25,000  ft. 
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Table  X 
of  Freedom 
Servo  Lock 


Listing  of  Computer  Output  from  a  Three  Degree 
Flight  Simulation  With  Adaptation  Loop  Using 


TIKE-SbC 


T I Kt- StC 
.500000 
. 525000 
.550000 
.575000 
-600000 
. 625000 
.650000 
.675000 
.700000 
.725000 
T IKE- St C 
.750000 
-.775000 
.000000 
. 625000 
. 650000 
.675000 
. 900000 
.925000 
.950000 
.975000 


Flight  Condition  7.’  M 
Adaptation  Servo  Gain 


1. 


1.8  ©  25,000  ft. 


i-.u-Sl.C 

I'M-  TA 

THl  I'AOuT 

DAMPR 

DA'-" PER  GAIN 

AC  1 UATC  S 

.000000 

.500000 

. 000000  ' 

. 700000 

.000000 

•  0  0  0  0  0  G 

. 023000 

. 496w86 

-.324431 

. 700  0  00 

.000000 

. 000 000 

. 050000 

. 484054 

-.  634  78m 

. 700000 

.OOOuOO 

. OOOo  Ou 

.0  7500C 

.  464490 

-.  926028 

. 700000 

.000000 

. OOOUOO 

. 100000 

.43  792  9 

-1. 193ol5 

. 700000 

.000000 

.0000 Oo 

. 125000 

. 405005 

-  1 . 43  4  34  2 

.700000 

.000000 

.000000 

. 150000 

.  3664  4  1 

-  1. 644  364 

. 700009 

. 000000 

.  00000’' 

. 175000 

.  32  30 36 

-  1. 82  122=1 

.700000 

.000000 

. OCOO 00 

.200000 

.  275644 

-  1 . 962866 

. 700000 

.00 00 00 

. 000 000 

.225000 

.225  167 

-?. 067863 

. 700000 

.000000 

. OOOuOO 

POs 


THtTA 


Till-  I  ADO  r 


DAMPR 


0 A  -1 P E i<  GA  I  N.AC  TUATOR  °0‘> 


.250000 

.  1  725  32 

-2.  1  354  06 

. 700000 

.000000 

.000000 

.275000 

.  1  ldobO 

-2.  165315 

. 700000 

.000000 

.  0  0  0  0  (  > 

.300000 

.064546 

-2.  1  58025 

.089002 

.000000 

. OOOUOO 

.325000 

. 0 1 1050 

-2.  1  i  45  76 

. 700000 

.0  152  75 

-.000452 

. 350000 

-.04  0  /  12 

-2.  054858 

. 7G000O 

.u  1  52  75 

—.0094 90 

.375000 

-.090410 

-1.91 7665 

. 127705 

.0  152  75 

-.02  1  140 

.400000 

156543 

-  1 .  766  75  4 

. 181400 

. 02 > > 8? 

-.028262 

.  425000 

1  78529 

-  1. 5° 64  59 

.700000 

.04254  7 

-.  0372  72 

. 450000 

-.215658 

-1.3786  75 

. 700000 

.0 42547 

-.0498,59 

.475000 

-.247287 

-1.  1  4  7o 54 

.205708 

.04294 7 

-  .  u5  74 4  1 

ThETA 
.  2  72925 
.292264 
. 505229 
.51 1 750 
.  6  11/  75 
.506190 
.294656 
.  278  >00 
.257  46 
.  935590 
THETA 
•.206468 
.  1  7  72  1  1 
.  146  >32 
.115119 
-.  063520 
-.052639 
. 022729 
. 0055 79 
.  03  1  7  71 
0554  1  7 


THC  TADu T 
-. 90  I  I  44 
- .  64  65  15 
-. 388959 
-.  15  5684 
.  11  54  10 
. 345684 
.557601 
. 744048 
902575 
031  7  56 
THE  1  A  DOT 
1.13 1944 
202945 
24  5.157 
2.6  1  9  75 
253045 
2  2  1  02  3 
167525 
1 . 092940 
.  998.985 
890160 


1. 


1. 

1. 

1. 

1. 

1 

.  1 


DAMPR 
.  199693 
.  126940 
.  129633 
.  13  7643 
.700000 
.700000 
. 700000 
. 700000 
. 700000 
. 700000 
DAMPR 
.700000 
. 700000 
.700000 
. 700000 
.266386 
.20 1 998 
. 700000 
. 700000 
.348036' 
394V69 


DAMPER  GAIN, AC IUAF0R  POS 


.084/35 
.067462 
.08  I  7  89 
.096048 
.110107 
.110  1  0  7 
.  1  I  0  I  0  7 
.110107 
.110107 
.  I  10 i 07 


- . 0565 64 
— .  052  ;>7>0 

9  .  04  7  I '  > 
-.039U44 
-.0263  I  1 

-.007646 
.Olio  '4 
.  0406  7  7 
.064  I  90 
. 0  8  49  2  ) 


DAMPER  GAIN  ACT UA I  Ok  POS 


.110107 
.110107 
.110107 
.110107  • 
.110107 
.  12094  7 
.  13  339 7  • 
.  1  33397 
’.  1  33  39  7 
.142196 


102290 
.116  15  7 
.  1266  9  I 
.  1  3  38  32 
.  I  3  7936 
.  I  39334 
.141761 
.  146842 
.  14  7  9  54 
.  142  796 


I  M  E  —  S 13  C 

THETA 

THE  TADOT 

DAMPR 

DAMPER  GAIN 

ACTUATOR  POS 

1 .000000 

.  076203 

.  771047 

.259582- 

.  149827 

.154114 

1.025000 

.0939 17 

.  6449  19 

.244  1  28 

.  16  084  3 

.  1238  15 

1.050000 

.  1084  19 

.514754 

.246414 

.  1  72239 

.  1 1280  7 

1 .075000 

. 1 1 9646 

. 383499 

,27-303  1 

.  183579 

.099 1 70 

1 .  1000  00 

.  12  7o  10 

.254555 

.28974 1 

. 194253 

.0837  15 

1 . 125000 

. 132406 

.  13  06  86 

. 700000 

.204510 

.0655  1  1 

1 .  150000 

.  1342  10 

.015725 

. 700000 

.204510 

.044834 

1  .  1  7  5  0  0  0 

.  1352  79 

-.  087  7  1.8 

.700000 

.204510 

.022451 

1  .200000 

. 129927 

-.  1 7  7  oOO 

. 700000 

.204510 

.000005 

1.225000 

. 124508 

-.•2  5  29  89 

. 700000 

.204510 

-.0208 69 

I  RE-SEC 

1  HET  A 

THE  1 ADOT 

DAMPR 

DAMPER  GAIN 

ACTUATOR  POS 

1 .250000 

. 1 1 7390 

3  135  70 

.  70 00 Oil 

.204510 

-.059 1 26 

1 .275000 

.  108/36 

-. 359950 

. 700000 

.204510 

-.05451  7 

1 . 300000 

.099493 

-.392683 

.700000 

.2045  10 

-.066542 

1. 32  5000 

08  9  SOb 

-. 4 1 3580 

. 700 000 

.204510 

-.0  75  S  70 

1 . 350000 

.078907 

-.  4  22  5  78 

. 700000 

.204  5  10 

-.08  16  15 

1.375000 

.,068832 

-.  it?  109  1 

.433248 

.2045  10 

-  .  0  8  5  5  S  1 

1 .400000 

.•05  789  1 

-.41 1/83 

.396940 

.211179 

—  . 0868  1  S 

1 .425000 

.  04  7  l  90 

-. 594597 

.220489 

.218755 

-  .  08  7  0  5  5 

1 . 4  50000 

. 0587  12 

-.  5  70  87  7  , 

. 185767 

.250743 

—  .  0ti6o  <  -! 

1.475000 

. 0795  16 

-.340152 

.  169  195 

.243599 

-.085652 
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Litis  r.iCiin  va'ue  which  h«.s  Scan  tabulate!. 


Table  XI  -  Comparison  of  Steady  State  Damper  Gain  Generated  by  the  Adapta¬ 
tion  System,  with  Theoretical  Gain  Required  for  -  0,7  at 
Each  Flight  Condition,  to  &  Sasic  Airframe  Damping 

Ratio  also  Tabulated. 


Adaptation  Servo  Gain  1.0 


Flight 

Condition 

Basic 

( 

Time  to 
Adapt 

Theoretical 

ICd 

Experimenta 1 

Kd 

1 

M  =  0,9 
10,000  ft 

0.440 

1.4  Sec. 

0.  156 

0.16 

2 

M  «  1.2 
10,000  ft 

0.177 

2.2 

0.303 

0.29 

a 

J 

M  =  1.4 
10,000  ft 

0.158 

2.1 

0.284 

0.28 

4 

M  =0.9 
25,000  ft 

0 . 244 

1.7 

0.2  49 

0,26 

5 

M  =  1,2 
25,000  ft: 

C.  100 

2.5 

0.372 

0.36 

5 

M  =  1.4 

25,000  ft 

0. 105 

2.5 

0.384 

0.38 

7 

M  =  1.8 
25,000  ft 

C.096 

2.5 

0.344 

0.34 

8 

Mb  0.9 
45,000  ft 

C.  141 

2 . 2 

C.  405 

0 . 44 

9 

M  =  1,2 
45,000  ft 

0.066 

?.  1 

0.538 

0.56 

10 

M  =  1.4 
45,000  ft 

0.057 

2.9 

0.503 

0.51 

11 

M=  1.3 
45,000  ft 

0.058 

2.8 

0,467 

0.47 
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It  is  seen  that  the  -by  taticn  times  coriare  l^Vurally  with  ^hcst  ol  the 
MIT  system,  and  as  would  be  expected*  these  flight  conditions  v. ith  1ower 
values  of  basic  damping  ratio  exhibited  longer  adaptation  times*  The  ex¬ 
perimental  values  of  steady  state  damper  gain  couM  be  determined  at  best 
with  two  place  accuracy,  but  even  so,  reasonable  agreement  with  the  theo¬ 
retical  values  was  obtained* 

Considerations  in  the  Practical  Mechanization  of  an 
Airborne  Digital  Damping  Ratio  Computer 

perhaps  the  best  method  of  reviewing  the  finally  "optimized”  Adapt¬ 
ation  process  is  to  follow  through  the  flov;  diagram  of  Figure  31*  v?hich 
illustrates  the  program  sequence  to  be  followed  by  the  proposed  airborne 
damping  ratio  computer.  The  primary  difference  between  the  airborne  com¬ 
puter  and  the  §-uJn  solution  programmed  into  the  three  degree  of  freedom 
flight  simulation  is  the  requirement  for  an  analog  to  digital  converter  for 
pitch  rate  input  signals  and  a  digital  to  analog  converter  for  output  signals 

to  the  adaptation  servo.  It  was  found  from  the  simulation  program  that  simj  le 

#  •  0 

trapezoidal  differentiation  of  0  and  0  yielded  sufficiently  accurate  values 

•  •  •  •  •  . 

of  0  and  G  for  use  in  the  calculations. 

From  a  short  preliminary  survey  of  computer  requirements  it  uas  deter¬ 
mined  that  a  special  purpose  digital  computer  having  a  130  word  storage 
capacity  and  a  repertoire  of  12  arithmetic  and  manipulating  instructions 
would  suffice  for  this  application.  And,  if  it  is  assumed  that  each  in¬ 
struction  requires  5  microseconds  for  execution  (a  representative  state- 
of-the-art-estimation),  the  sequence  from  raw  data  input  to  generation  of 
a  new  adaptation- servo  command  will  require  12*600  microseconds,  or  in 
other  words,  the  total  cycle  will  have  a  maximum  possible  repetition  rate 
of  83  times  per  second*  Hox;ever,  the  cycle  utilized  in  this  study  was  re¬ 
peated  only  40  times/second,  which  leaves  a  comfortab1©  execution*  s^eed 
margin. 


8? 


Figure  31.  Flow  Diagram  of  the  Proposed  Damping 
Ratio  Computer  for  Use  in  the  Adaptation  Loop 
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VII  i  COifcLUSlOitt 


It  has  been  demonstrated  that  a  simple,  special -purpose  digital  com¬ 
puter  can  be  utiMzed  to  determine  the  natural  frequency  and  damping  ratio 
of  the  short-period  mode  of  the  free  longitudinal  oscillations  exhibited 
by  a  representative  high  performance  fighter  aircraft.  And*  it  has  also 
been  shown  that  the  calculated  values  of  damping  ratio  can  be  used  in  a 
closed  loop  manner  to  vary  the  gain  of  a  self  adaptive  pitch  damper  which 
is  capable  of  maintaining  constant  airframe  short-period  damping  ratio  over 
the  full  speed  and  altitude  range  of  the  vehicle.  The  system  investigated 
is  capable  of  performing  the  adaptation  process,  from  basic  unaugmented 
airplane  to  £  =  0.7,  in  a  period  of  time  ranging  from  1,4  to  3.1  seconds 
depending  on  flight  condition;  and  this  speed  of  operation  was  shown  to  be 
similar  to  that  of  the  most  recent  MIT  flight-vehicle  self-adaptive  system. 

In  conducting  an  evaluation  program  with  the  proposed  mechanization, 
in  conjunction  with  a  longitudinal  three  degree  of  freedom  flight  simula¬ 
tion,  it  was  found  that  measurements  of  pitch  rate  rather  than  pitch  angle 
were  necessary  to  extract  the  short  period  characteristics  from  the  combin¬ 
ed  phugoid  and  short  period  motions.  First  and  second  derivatives  of  pitch 

rate  are  then  determined  in  the  adaptation  computer  and  sampled  values  of 

*  *  •  •  • » 

0,  9,  and  0  set  into  the  iterative  routine  for  determination  of  short  period 
damping  ratio  and  natural  frequency.  The  system  requires  no  perturbation 
signals  or  limit  cycling  for  operation  of  the  identification  process*  how¬ 
ever,  the  system  will  adapt  only  while  the  airframe  exhibits  pitching  motion. 
It  should  be  emphasized  that  the  mechanization  investigated  in  this  study 
is  not  a  complete  longitudinal  control  system,  in  that  there  is  no  pro¬ 
vision  for  external  pilot-input  signals.  However,  it  was  proposed  in 
Chapter  V  that  Corbin* s  forward- loop-gain  computing  system  could  be  con¬ 
verted  from  analog  to  digital  form*  and  used  in  conjunction  with  the  present 
computer  routine  to  enable  the  introduction  cf  external  .ontrol  inputs. 
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- 1 hv  1  1  ir  a  study  of  1  e  •  .e  0  forcing 


the  rcmviyr^euce  of  *  <*t  sextos  of  £-0^  Iterative  calculations  It  is  felt 
chat,  if  the  calculations  could  be  made  100  percent  converger tv  the 
speed  arid  accuracy  cf  the  adaptation  process  could  be  greatly  enhanced » 
Other  factors  which  vould  require  investigation  before  a  practical  mechani¬ 
sation  could  be  accomplished  are:  1.  the  effect  of  noise  (atmospheric  tur¬ 
bulence)  introduced  into  the  forward  loop,  2.  the  effects  of  control  sur¬ 
face  rate  end  amplitude  Uniting^  and  3*  the  influence  of  transducer  in¬ 
accuracies  and  noise  on  the  operation  of  the  adaptation  loop. 
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+  .  &&  O 

Cl a 

4-.4Q 

3.  eo 

3.o9r 

2.. 2.9 

Cut 

.  008 

.  3SQ 

.  ZG  + 

.  /  70 

Cr>e 

.2-34- 

.  2.4-1 

.  )S8 

-  0*5 

Cn6r 

3  70 

-.3  2-8 

2S6 

-.13  8 

Cngct 

.  0131 

.  ooze 

.002.8 

.  OO  ZB 

Crip 

4-72. 

+08 

-.4-7  0 

-.  a  7^> 

C  np 

•  too 

.  075 

.  O  70 

.0  62. 

Cj.0 

-i+s 

-.75 

-.13  0 

037 

Cl-Sr 

.  1  Z  + 

.  // 5 

.too 

.  O  7*5“ 

.on  6 

.  OOSJG 

.  O04OZ 

.OOJ36 

CjLy. 

.  163 

./S’O 

.738 

.114- 

CjLp 

3  BS 

-  3>/ 

-.2.6 

-  a 

Cye 

-1.  ^  2. 

/ 

-1.0 

8S  '  ~ 

c*. 

- 

Cy  Sr 

.  80  & 

.  soo 

'  ▼ 

.3^r<7 

Cl, 

CDo 

.  • 

,  t 

f 
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*> 


Stocbilit 


7 


Deri  vat) 


/ ves 


Altitude  fS.OOO  ft. 
tA~  O.Q  J.Z  A1- /.4-  M-  /•  8 


cM;t 

-.39/ 

-.  74-5 

-  .555 

-  .367 

Cm  y 

-4-.de 

-4.  7/ 

-3.83 

-r-  2.90 

Ci”  ex. 

-.72.2. 

-  /.  JG 

-  .9/5 

-  .GO 2 

ClTl^ 

-  A3  G 

+  2.37 

4-1.0  7 

■f  .  90 

Cl  a. 

4-.  74- 

4. 19 

3.43 

3.4-7 

CcU 

.  725 

.47/ 

.352. 

.2-32- 

Cue 

.  2G7 

•  Z07 

.  zoo 

.095 

CnSr 

-.399 

-.  364 

-.3  61 

~.  2.0/ 

Cnja. 

.0/3/ 

.  0/0 

.o/OZ 

. 0/0 

Cnr 

-  .54  o 

55/ 

- 470 

32.3 

Cnp 

.12. 

.  10  2. 

.0  84- 

.082. 

Cjte 

-.156 

-.17/ 

—  1  60 

-.0  95 

Cj5r 

•  /  34- 

.  1  ZB 

.  J13 

.  OQ7 

Csia 

.  0/33 

.0127 

.0008 

.00  9 

CjL  p 

'.  2.4-2. 

.176 

.136 

.114- 

Cjl  p 

-.42 

-.36 

-  3/5 

Cya 

3  o 

- 1.J3 

- 1.10 

-  .5  7  ' 

Cyr 

- 

Cy  £  r 

.  8  70 

.  7/0 

•  56  8 

.4/9 

CyP 

Co* 

.  * 

77 


( 
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Table  cl  Airplane  Physical 


c ha ra c t er  1st  1c ? 


Airplane  macs  1,110  slugs 


Tx 

19,542  slug- ft. 

ty 

129,  5.33  slug  ft,2 

lz 

142,  398  slug- ft } 

Txz 

8,219  slug- ft. 2 

Wing 

area 

394.04  ftd' 

Wing 

span 

96.53  ft. 

Wing 

mean  aerodynamic  chord 

12.24  ft. 

Wing 

aspect  ratio 

3.39 

I  , 

;{e 

13,93  slug- ft. 2 

a)e 

1*259  rad/sec*  (Cruising) 

110 


A  S?rmn  Cl-  f'K  TCI.  oi’AL  L!.Ti lLV  A  r\’£i 


DERIVATIVE 

CGMElClS 

C16  a 

Rolling  moment/  aileron  deflection  lit",  le  import 

in  dynamic  stability* 

Clfr 

Rolling  moment/  rudder  deflection  cross  coupling  fcerr. 

of  little  importance  In  dynamic  stability. 

C  7 

f 

F.c  1 1  log  moment/  rol  ■  rate  -  viscous  damp  trg  In  rol1, 
balances  aileron  moment  to  produce  r,  constant  lolling 
velocity  for  a  given  aileron  def ! 'ction*  has  small 
damping  effect  on  lateral-directi onal  spira1  mode, 

CI 

r 

Rolling  moment /  yanking  velocity  -  contributes  slightly 
to  Dutch  rol1,  most  important  cause  of  spiral  divergence 
mode* 

s 

Rollirng  moment/  sideslip  angle  -  dehedral  effect,  large 
negative  value  v/ill  produce  good  spiral  characteristics 
and  poor  Dutch  roll  characteristics . 

°l0 

Rolling  moment/  sideslip  angular  rate*  usually  negligible 
in  magnitude. 

C 

in. 

Xt 

Pitching  -Moment/  elevator  deflection*  control  derivative, 
little  importance  in  dynamic  stability. 

C 

ra 

q 

Pitching  moment/  pitch  rate,  viscous  damping  in  pitch* 
large  effect  on  damping  of  longitudinal  short  period 
mode. 

c 

*05 

Pitching  moment/  rate  of  change  c£  angle  of  attack*  due 
to  aerodynamic  time  lag  and  aeroelastie  effects ,  contributes 
to  damping  of  longitudinal  short  period  mode. 

c 

mac 

Pitching  moment/  angle  of  attack  -  produces  lestoring  moment, 
establishes  frequency  cf  longitudinal  short  period  mode. 

C 

rr.  p 

Pitching  moment'7  sideslip  angle  -  aerodynamic  cross  coupling 
derivative,  usually  negligible  in  magnitude* 

C 

n<5  r 

Yawing  mos^nt /  rudder  deflection  «  control  derivative-*  little 
importance  in  dynamic  stability* 

C 

n 

r 

Yawing  moment/  yaw  rate  -  viscous  damping  in  yAw^  major  contri 
butcr  to  damping  In  Dutch  roll  (lateral -direct  tonal  short  peri 
cd  mode^ „ 

111 

n 

r 


■VzLlVi.t  !  V Ci  4 


r  1 1  tv  t  -  **  .-.ro^ynauiic  cross  ecupHng 
cnt r  routes  slightly  to  damping  In  Ditch  roll, 


C 


v 


fax.Lg  .>  ant/  sidesKi  angle  -  produces  directional  re¬ 
storing  T:?«.-.3nt  ,  establishes  frequency  cf  Dutch  roll  node* 


C  . 
n0 


Yawing  moment/  sideslip  angu1ar  rate  -  contributes  slightly 
to  damping  cf  Dutch  roll. 


C 


Yawing  r..c*ment/  aileron  deflection  -  control  cross  coupling 
derivative,  little  importance  in  dynamic  stability. 


C.r  cide  force/  sideslip  angle  -  contributes  to  damping  of  Dutch 

‘fi  rolU 


Gv  Side  force/  roll  rate  -  aerodynamic  cross  coupling  deriva- 

~p  tive,  magnitude  usually  negligible® 

Cv  Side  force/  yaw  rate  *  magnitude  usually  negligible® 

r 


Cv  .  Side  force/  sideslip  angular  rate  *  magnitude  usually 

A  &  rtegl  igib  le . 


r> 

Li 


cc 


c 


'CL 


c 


Lift/  angle  of  attach.  -  effects  frequency  of  longitudinal 
phugoid  :iode,  ha^  wi jor  effects  on  aircraft  maneuverability. 

Lift/  rate  of  change  of  angle  of  attach  -  due  to  aerodynamic 
time  lag  and  agroelastic  effects,  has  negligible  effect  on 
dynamics „ 


Lift./  e1evator  deflection  *  magnitude  negligible 
ticnal  aircraft,  has  effect  on  tailless  aircraft 
with  short  tail  length. 


for  convert- 
and  those 


C 


D 


Airplane  drag  coefficient  -  primary  cortributsr  to  damping  of 
longitudinal  phugoid  oscillation. 


Drag/  angle  of  attach  -  negligible  effect  or* longitudinal 
phugoid  mode  because  this  mode  takes  place  at  constant  angle 
of  attack. 


A*  1  hN  *1  a  IV 


FQ  Allijf.i  Cr  l*OtIO;. 

The  equations  r  ?  c o  are  leif.ed  from  th*j  cvo  far  II  i  -  press- 
tons  :  Force  ~  '  x  acceleration  and 

Moment  =  Wou.ent  ■:  f  Inertia  x  Angular  Acceleration 
From  classical  mechanics,  the  equations  of  motion  cf  a  rigid  body  in  free 
space  with  no  aerodynonic  or  gravity  forces  acting  >re 

/*;(  =  m  (  u  f  -  rv) 

/-y  -  m  (  v  +  ru  -  pw) 

-  m(w  t pv  —  u ) 

Z=  +<^r(Ij£  ~Jy)  ~  +**)  +~£)(y  (rp~<%)  +  Iy^(rZm^cF) 

M-  ^  Iy  +  (p$  -n)lyi  ~  (p  -f-yr)  lxy  -h  rp  (Iy  -Iy.)  +  (p^-r2)  IXi 
N=  rli  i-  (yr  -p) IXi  -  (y-hrp)Iyz  + pc^  ( Iy- Tx)  +  (cf-p)  Ixy 

Note  that  all  symbols  ere  def?iiei  in  figure  1,  It  is  seen  rK:c  the  ’ast 
two  terms  in  each  of  the  first  three  equations  are  th:  pi  civet  of  an  angu¬ 
lar  and  a  linear  velocity  and  thus  represent  cenctipeta1  c tele rations  dir¬ 
ected  along  che  respective  ^xes.  The  1st  four  t«ros  f.n  e  •:  l  the  last 
three  equations  represent  nc^ents  created  by  inertial  cross  couplings  These 
terms  arise  fren  a.  ;v:.$  tty  of  sources  including  gyroscope.  •  rt-  es^ior.o  1  torques 


e filial  couH -.i:g  1  1  1  _  g?.v&».  in  *  •i'?  $  •-  £ 

pointed  cut  tftifc  d*c<$  cer  r  \x)  r>oi  stfrses  l^ccv -  n  tv  *s*strous 

results,  when  the  a-  •  Mtulss  .  ivl  vales  of  aircraft  »--»otirn  1  a  t  ge  (ref. 

3).  It  is  recalls  l  that  tbn  a  It  era  ft  possesses  vert  ( 1  ditch 

causes  the  products  of  inertia  Ivv  an  1  I./v  to  oe  ter..  ...  This  e  I  .  a  tnaces 
two  inertia  ^  coupling  terms  cron  each  equation,  IVext*  the  coim  ottnes  or 
force  due  to  gravity,  and  the  qcro  iyna-ilc  forces  and  on  eots  are  included. 


P  5  _ 

ty) u  —  rv  -  w  —  /ry  *s/7?  0  ~  V^vS  Undine  Thrust 


mv 


=  yc ?w  - ru  nhQ  cos  O  s'm0  +  Cy  VL & 


p 

m\N  -  <^u  -jov  "  mcj  cos&  cos(p  ~  j"  '/ ■“ £ 


plx  -  yr  (Iy-I?)  +  rl^  +  py_  +Cj-gV-66 


yly-  pr(ll-ly,)  -p^Ixt  ~  Iy.*>  ^ a  r  ^  Cm  2.  VLSc 


lXi'P  lX  't  ~  lXe 


p  ^ 

p<^  (l\~ly) +■  p  Ixi  ~  %r  ?xi  i~^Xe°^^%  ~jT  V  Ch 


The  I  CP  terms  rar  reset,  t  he  snjuV.r  nortneutn  of  rotating  wires, 

Xa  a 

which  when  Haiti*  lied  by  a  yaw  or  pitch  singular  -8? ocU>*  j  *  gyro- 

sccpic  process  lorn  1  torque  about  the  pitch  cr  yaw  ax^s  rtcjc-ct"\ely*  These 
terms  are  significant  in  a  turbojet  powered  aircraft.  It  should  be  noted 
that  the  wingspan,  b-  is  u ce*1  nr  v*c  -ret  t  &v:  to-.  th«  •Krrbw;.  ^  .  maw ts  io 
the  roll  and  yaw  equations  whereas  .ear  aerodyne  ib'c  d-oc  T„  *  -  used  as 


iv. 


f  {'  \  ..  *r  ici'  _*•  i  •:>. 

I  c  .*ui >  X.  ?  .  '!iw  -u/  ulo,r 

ecr’^onencs  *'«f  1.5.. -t*:  *  =  *  u,  v  .  'av*  > 

into  nova*  "a  .■  •  1  i  &  t  4  •  qi- '  4 1  ies  „ 

U  ~  V  }  oLirjoloLn e  Forward  velocity 

\/V  V  6in  3  **  VB 

V  Sin  (X  V  CL 
v  ^  \/f? 

Vac 

Where  a:  and  &  ,  the  aerodynaff*ic  attach  Juv’  siAifc?  ift 

inclination  of  th«  forvar  •  velocity  vector  to  the.  T  bci'v  $x;.s*  The 
complete  equations  of  motion  with  the  aerodynamic coefficients  T.  Expanded 
form  are  shown  in  table  I,  The  reader  will  see  that  the  aerodynamic  ro¬ 
tary  derivative  terms  are  multiplied  by  certain  add \t tonal  flight  Condition 
parameters  in  order  th*r  the  value  of  the  derivative  cmv  term  in  as  constant 
as  possible  over  a  wile  range  of  flight  conditions*  example  1°  uae 

of  the  term  with  derivatives  taken  with  respect  to  roll 

terms  for  other  derivatives  are  listed  in  the  table  of  sy^hoSs* 

One  final  set  of  equations  must  be  used  to  enable  the  conversion  of 
angular  rates  about  the  body  axes  to  angular  rates  with  caffpsct  to  the 
earth  coordinate  system •  The  necessary  Suler  equations  are  iref.  5. ) • 
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0  =  ^  cos  cp  —  r  sin  (fi 

<p  -  p  r  c?  cin<p  tan  <p  -h  r  cos  <f>  tan  Q 

¥  =  ( %  sJn  <f>  +  r  co 5  $ )  sec,  O 

•  • 

G  and  ^  can  be  integrated  to  obtain  0  and  0  which  are  used  in  the  equa¬ 
tions  of  mot  ton «, 
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DSFJVATION  OF  P‘E  LO NO  I T  i3  i  i  I  *  A  L  SHORT-VERTor  iuu!S  Ffik  FVIieTlG.;  .'0  IHK 

EQUATIONS  OF  MOTION 


Of  the  two  longitudinal  modes  of  airplane  motion,  the  phugoid  and 
the  short-period,  only  the  latter  mode  is  of  interest  in  stability  augmenta¬ 
tion  considerations  because  the  phugoid  normally  does  not  cause  any  deter¬ 
ioration  of  handling  qualities.  In  Section  Til,  a  transfer  function  of  the 
mode  of  interest  in  stability  augmentation  systems  must  be  included  in  the 
block  diagram  of  the  proposed  [itch  danger  to  enable  an  analysis  of  the 
system.  Thus,  the  longitudinal  short-period  transfer  function  must  be  de¬ 
rived  from  the  airplane  equations  of  motion.  fhe  Jegrees  of  freedom  of 
interest  being  pitch  and  normal  translation  for  which  the  corresponding 
equations  are  shown  below. 


>\  Ac  .  ,  AC  z*  r 

C7  =  ~  C  t£  ~h  — Ct 


Ac 


y 


£  yjy  ^  T_  r.  j  ^rr)^  Act 


'?**  Q  +  ly 


OC,  -  G  +  y  cos  6  C'Lcc  ^ 


The  pitch  damping  derivatives,  Cm^  and  ,  are  usually  added  together 

to  form  an  equivalent  Cm a+a.  >  engine  torque  and  tail  lift  terms 

have  been  neglected  because  of  their  lor  longitudes s  As  sum-,  t  ion*  neces¬ 
sary  for  the  derivation  of  the  transfer  function  are:  1.  s~  equilibrium 
condition  exists  a’or.g  the  *  axis  when  airplane  weight  »  a i x r. ' une  lift, 

W=  — '"’^'”0.(2,  <^-r^r"r)  *  ax)('  ('°’S  &  ^  l  • 

The  difference  between  CC  and  CC  -trim  will  be  denoted  Cc, 


1 17 


equations  of  wot  Lon  then  he  cor  e  : 


Ac 

iy 


Cnr,  't  +• 


A  G± 

2-VIy 


Cm%+a.  Q 


-h 


Ac 

Ly 


Cmcc  ^7 


01  =  0- 


A 


7nV  ^ L 


C,  _  a;. 


CO 


The  next  step  is  the  introduction  of  the  differential  operator*  S,  and  the 

0 

elimination  of  angle  of  attack  from  the  equations  in  order  that  the  — — 
equation  may  be  obtained. 


sLe  = 


Ac 

h 


Cm'ij.'t  + 


Ac± 

2  Vis/ 


Cm^ ri  06  +• 


Ac_ 

h 


Cm  a,  ax., 


ccb  -  OS  - 


mv 


Cl<x 


or 


96 


5  +■ 


'Cl 


CO 


Substituting  CC, 


os2-- 


Ac2-/-.  css  Ac 


2.  VI 


Cmy+a,  60 


ly 


r  66 

r._A_ 


6  KOV 


A  c  n  : 

Zy  C/77'i  £ 


0 

The  -j —  transfer  function  is; 

/f 


0 


Ac  / 

77  Cmi* 


my 


Clcc) 


!  t  c  /  Cu  x  C-A.C - A  c  —  AC2-  ~  a  S'  Ac  /-  \ 

T^mvci<z  zviy  ^.viy  nov  ('Lcc~~  ~^n^ccJ 


Zy 


Rearranging 


e  


AC 

/v 


Cmtt  (S 


^  r  ^ 

^LocJ 


mv 


'•t  '  s^+sfacL,-mcn%)-ly&-c„%±7cLa-&-c„I)) 


11.8 


‘:?h^  fi 


'  n  ■ 


>r  {;  __<<v  :  s  «-hn*  n  he  'V 


=  /recs^r*  j _ 

;£  S(SZ  i~2£uJr}5  i~ CUn *) 


fi$  I  s  the  a  Lj.’i  1  -r.ne  a1  evatur  e  f  feet  tveness 
7a  'i s  the  H  lyht  path  t*  ime  constant 

£  Is  the  short- rei'iod  damping  rot  i  o 


(j)n  Is  the  short-period  natural  frequency 
It  can  be  seen  that  the  steady  state  transfer  function  is  essentially 
a  pitch  rate  command  system  for  a  given  horizontal  siaMI'ger  deflection* 


IQ 


NOT  7  CI\T  T^S 


:IC  SYSTEM 


A1  though  analog  cot*ute:es  ere  general  ’y  used  to  simulate  system 
dynamics,  the  systems  investigated  In  this  study  v-ere  simulated  ou  a  stored 
program,  general  purpose  digital  computer  using  techniques  similar  to  ana¬ 
log  computer  programs.  The  techniques  are  similar  in  that  both  methods  use 
some  form  of  integration  combined  v?ith  other  operations  such  as  summation, 
subtraction,  and  multiplication  to  solve  a  set  of  simultaneous  differential 
equations  of  motion;  each  equation  of  motion  representing  a  degree  of  free¬ 
dom  of  the  system.  An  analog  computer  program  may  he  regarded  as  a  continu¬ 
ous  parallel  solution  of  the  equations  of  motion,  whereas  the  digital  method 
employs  a  step  by  step  serial  solution,  A  number  of  digital  integration 
techniques  are  available,  ranging  Cron*  simple  rectangular  integration  to 
highly  accurate,  but  time  consuming  methods,  such  as  the  Runge-Kutta  numeri¬ 
cal  Integration.  The  technique  selected  for  this  study  is  a  modified  form 
of  the  trapezoidal  rule  and  is  shown  below, 

X'-,  +  j-E  3  Xn_, 

vhere;  n-h  denotes  values  computed  in  previous  cycles 

h  is  the  width,  in  time,  of  the  computation  interval, 

A  major  difference  between  the  digital  solution  and  the  concept  of  the 
analog  flow  diagram  lies  in  the  fact  that  the  digital  method  perforins  in¬ 
tegration  of  lower  ordered  derivatives  first.  As  an  example,  the  steps 
in  the  computation  cycle  of  a  seconl  order  system  solution  are  shown  he^ow: 

X/?  =  X/?-/  ^  ^X/7-/  ~ 

X/7  -  X/7-/  +  ESXv-j  ~  X/) 7 

Xn  =  -d^jJn-CUf  Xn 
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The  cot.  j*u tat  Ion  eye1*.  is  repeated  at  a  high  rate  and  the  values  of  X^jSX 
generated  In  each  pass  rep  resent  joints  on  a  time  history  of  system  motion* 
T7hen  a  more  complex  system  having  a  greater  number  of  equations  of  motions 
Is  simulated,  the  integrations  for  each  equation  of  motion  arc  performed 
in  a  serial  manner  such  that:  one  pass  is  executed  through  all  equations 
before  recycl ing. 

The  accuracy  of  this  technique  is,  of  course,  a  function  of  the 
vidth  of  time  increment  between  computing  cycles,  and  it  was  noted  that 
when  using  the  method  described  above,  integrator  errors  were  manifested 
as  an  undamping  of  the  system  dynamics.  In  some  lightly  damped  cases  it 
was  found  that  cumulative  integrator  errors  actually  caused  the  solution  to 
diverge,  the  corrective  action  being  the  use  of  finer  calculation  inter¬ 
vals.  It  has  also  bean  noted  that  the  required  fineness  of  calculation 
intervals  is  a  function  of  both  the  period  of  the  system  mode  of  highest 
frequency  and  the  damping  of  this  mode.  A  rough  rule  of  thumb  is  the  use 
of  200  divisions  per  period  for  $  =  0.1,  down  to  50  divisions  per  period  for 
£  =  .5,  to  obtain  one  percent  solution  accuracy. 

The  following  five  pages  are  sample  listings  of  two  of  the  FORTRAN- 
programming- language  programs  written  for  this  study.  The  first  is  a 
second  order  system  and  the  second  a  three  degree  of  freedom  flight  simula¬ 
tion.  The  programs  were  run  on  both  the  GDC  1604  and  the  IBM  7090  general 
purpose  digital  computers. 
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SECOND  ORDER  PLANT  IDENTIFICATION  ROUTINE  3 


M=i 

1  -READ-  ILiPUT  TAPE  3,2,  ACHM, 
1A,STHETA,  ALPHAT,  THRUST 


ALT,  DENS  TV ,  VELX,  DYNPRS,  TREJA,  C.THfcT 


1  ,  -5  » 


:  FORMAT  (F4.I,  16, 

READ  INPUT  TAPE  3 
J- FORMAT  (  6F6.3,  F8.4) 
WRITE  OUTPUT  TAPE  2,4 

t — EO  RM  A  T _ C.  ITiiy  /  5  8H  _ 

1NUKB  5R  =  F4.i,  I2H 


F8.7,  F7.2,  F8 
CMIT,  CMOA, 


1  ,  4F7.6,  F8.2) 
CMA,  CLA,  CLDLLE, 


CDU,  TIIEDD 


INITIAL 
J  =  1 
Ji  =  l 
K  =  1 


CONDITIONS 


,h,  ACHM,  ALT 
-THREE.  DEGREE  OE-  .FREEDOM 
ALTITUDE=  16) 


FLIGHT  SIMULATION  .MACH 


_ L  =  1 _ 

L  1  =1 

-  T IME=G • 

SR  VGN= .25 
_DAMPR= . 7 
DMPGN=0, 

- DAMP.R2  =  .4 _ 

ACCLN=0. 

OM£GAN=0. 

50  A L PHA= . 5 

_ ALPHD=0.  _ 

AL  PHD!  =  0 . 

_ THETD_=0. _ 

THETDi =0 . 

- TH£DDi  =  THEDD 

THTTD=0. 
THTTD2=0. 
-VELCTY  =  VELX 

- V  ELXD=.Q.  _! _ 

VELXDi =0. 

_ A  C  TP  0  S  =  0 . _ 

ACTVL=0. 

_ A C TV L 1=0. 

AC  TC  L  =  0 . 

_ ACTCLI=0. _ 

S I GNAL  =  0 . 

AL  PHA2.=  *  5 
A  L  PHD2  =  0 . 

- _ THET A2  =  .  5 

THETD2=0 . 


_ T H  ED  D  2  =  THE  DD 

5  WRITE  OUTPUT  TAPE  2 ’,6 

—6  FORMAT  ( 105H  TIME-SEC  THETA  THETADOT  FWD  V 

1UAT0R  POS  NORM  ACCEL  DAMPER  GAIN  OMEGAN  J  DAMPR) 

7  WRITE  OUTPUT  TAPE  2, 8, TIME,  THETA,  THE TD , VEL CTY ,  ACTPGS 
1DMPGN,  OMEGAN,  DAMPR 

8  FORMAT  (  9  F  1 2  . '6  ) 

I  F  (L-1000  )  9,  12, » 12 

9  L  =  L  +  1 

IF (Ll-200) 51, 52, 52 

51  L 1=L 1+ 1  1 

GO  TO  14 
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ELOCITY  ACT 
,  ACCLN , 


SECOND  GRDcR  PLANT  I  DENT  I  F I  C  A1  I  ON  ROUTINE  3 


:  52  L  1  =1 

_  go  _m_5o..  - _ _ _  _ _ _  _ _  .  _ _ 

12  CALL  EXIT 

ACTUATOR  SIMULATION 

14  ACrPOS=ACTPOS+.OC23*(3.*ACTVL-ACTVLi  ) 

ACTVL l=AC  T VL 

ACTVL=ACTVL+.Q025*(3.*ACTCL-ACTCLI ) 

_ ACTC  Li  =  AC-T.CI _  .  _  _ 

ACTCL=-ACTVL*49.-ACT  POS* 122 5.+ SIGNAL *122 5. 

INTEGRATION  OF  LONITUDINAL  ACCELERATION  TO  OBTAIN  VELOCITY 
VELX=VELX+.0025* ( 3 . *VELXD-VELXD1 ) 

_ VELXD1=VELXD 

VtLXD= (THRUST/1110. )- ( AOVtRM*CDO )- ( AQVERM*( (CL A* ALPHA) **2.  )/(  3. 14. 

_ 15 9 *3 ..3.9)  )rJ32L._18*SiNF(  THETA )  )  . . 

Vc  LC  TY  =  VE  L  X 

.  ..CALCULATION  OF  AERODYNAMIC  PARAMETERS  DEPcNDENT  ON  VELOCITY 
DYN.PRS  =  DENSTY*(V£LCTY**2.  )  *394.34 

_ ACIYY=DYNPRS*12. 24/129533. 

AC2VIY=(ACIYY*12.24)/(2.*VLLCTY) 

_ AMV=D YNPR S/ ( VELC TY*1 1 10.  )  _ 

AOVERM=DYN PR S/1110. 

SOLUTION  OF  EQUATIONS  OF  MOTION 
THETA=THETA+.0025* (  3  .  * TH ETD-TH ET D1 ) 

'  THFTD1 =THFTO 

THETD=THETD+.0025*(3.*THEDD-THEDDI ) 

_ THEDD1  =  THEDD. _  ...  _  _ _ 

THEDD= ( AC  I YY*CMI T*ACTPOS )  +  ( AC2 V I Y *CMQA*TH ETD )  +  ( AC  I Y Y *CMA* ( ALPHA-AL 
_ 1PHAT ) ) 

ALPHA= ALPHA*. 0025* (3. *ALPHD-AL PHD! ) 

ALPHOl=ALPHD 

..  ALPHD  =  THETD+  (  32. 1 8*C0SF ( THETA ) /VELCTY )- (  AMV*CLA*ALPHA) 

_ CALCULATE  NORMAL  ACCELERATION  _  .  -  _ li.  _ _ 

ACCLN=(CLA*ALPHA*DYNPRS)/ ( 1 1 1 0 . *32 . 1 8 ) 

_ PITCH  DAMPER  SIMULATION  _ _  _._V 

S I GNAL  =  THE  T0*DMPGN  .  .  •  •  - 

_ lTIME  =  TIME  +  .005  „_v> 

DIFFERENTIATE  PITCHING  ACCEL'ERAT  I  ON 

_ THTTD= (THtDD-THEODl ) /.005 

ADAPTATION  SERVO 

DMPGM=DMPGN+ ( ( . 7-DAMPR ) *SRVGN* .005 ) 

1  I F  ( J-5  )  15, 16,16  •  ' 

15  J=J+1 

GO  TO  14 

1 6  _  J  =  1 

CALCULATE  NATURAL  FREQUENCY  AND  DAMPING  RATIO-  ASSUMED  DAMPR=.4 
DAMPRi  =DAMPR2  ‘v 

21  GAMMA=DAMPR1*THEDD2/THETD2 

BtTA*SQRTF (GAMMA**2.rTHTTD2/THETD2 ) 

OM EG  AN =-  G  A  MM  A  + 1>  E  T  A 

DAMPRI =-  (  THT'l  D/(  2.  * OKtGAN* TH EDO )  )  -(  (  OMEG AN*TH2TD  ) / ( 2 . « THEDD) ) 

27  IF(J1-20)30»3‘1  ,31 

30  J1=J1+1 
GO  TO  2i 

31  J 1 = 1 

I F ( OMCGAN-8. ) 25,26,26 
25  IF  (  0  M  E  G  A  N ) 2  6 , 3  2 , 3  2 
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l 

32  I  F  (DAMP R 1-1.  )  33,26,26 

,33_IF  (.DAMP.R1 )  26  ,.3-4.,  34  _  . . . .  . . 

34  DAMPR  =  DAMPR.l  "• 

DAMPR2=DAMPR1 
GO  TO  29 

26  DAMPR1 =DAMPR2 
3  5  GAMMA  =  DAMPRi*THEDD?>/THETD2 

.  BETA=SQRTF  { GAMMA**! . -THTTD2/ LHETD2  )  _ 

OMEGAN=-GAMMA-oETA 

DAMPR1=- ( THTJ0/(2.*0MtGAN*THEDD) ) - t ( OMEGAN*THETD)/{2.*THEDD))_ 

36  I.F(  J1-20)3T,  38,38 

37  J 1 =Ji+ 1  .  “  -  . .  _  __ 

GO  TO  35  % 

„38„Jl  =  i  .  .  _ . _  _  ___ _  _ 

IF (OMEGAN-8. ) 39,40,40 
39  IF(0MbGAN)40,4l,41 
41  I-F  (t)AMPRl-i.  )  42 ,40,40 
.42  IF  (DAMPR1  )  40,34,34 
AO  DAMPR= . 7 

, _ wR.iTfc_QUTPUT_  T APE  2,43 _ _ 

’43  FORMAT  (31H  NO  SOLUTION  FOR  DAMPING  RATIO) 

29  THETD2=THETD  .  . 

THEDD2=THED0 

THTTD2=THTTD  .  .  _ _ 

IF (K-iO) 17,18,18 

GO  TO  7 

18  K  =  1  _ ■ _ _ _ _ - _ 

GO  TO  5 

E\D(1, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0 ) _ 
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\ 


XEG 

1  =  ] 

J=1 
K  =  1 

XB0=0* 

XCDO=C 
X  C  =  1  » 

XDN=0 • 

X  N  =  1 « _ 

T  I  ME=Q  • 

xcco=o. 

XCDOO=0* 

WRITE  OUTPUT  TAPE  2,1 
FORMAT  ( 1  HI //44H  PLANT 
C  AMPR=  „  4 

XDDN=-2*  4*XDN-36o*XN 


IDENTIFICATION  ROUTINE  -  SECOND  ORDER) 


XDINCR={ XODN+XCDO) **00025 
XCN=XCO+XDINCR 
XCDO=XCDN 

XINCR= (XDO+XDN ) *.00025 
XN=XN+X INCR 
XDO=XDN 

TIME=TIME+.0005 
IF! 1-1000)4,5, 5 

4  1=1+1 
GO  TO  2 

5  1  =  1  _  _ 

WRITE  OUTPUT  TAPE  2,6 

6  FORMAT  .(  1  H  //55H  _ _ TIME 

1EMENT) 

7“fcRMT0UhLI)PE  2,?’  TIME’-M0N’  - ' 

ITERATIVE  CALCULATION  OF  N AXUR.A  1 — FREQUENCY  AND  OAMPlNG  RATIO  OTVF 

A  STARTING  VALUE  OF  DAMPING  RATIO  =  0.4 - ~  RAT_1 0_  GIVE 

WRITE  OUTPUT  TAPE  2', 8 

8  FORMAT  ( 25H  GMEGA-N 

9  ALPHA=CAMPR*XDOO/XO 


ACC E  L  ER ATI ON V  E  L  OC I T  Y D I S  P  L  A 


DAMPRAT ) 


*XDN) ) 


B ETA= S CRT F{ ALPHA **2*-XDDOO/XO) 

OMEGAN=-ALPHA+BETA 

DAMPR  =  -(XDDN/ (2*  »OMFGAN*XDN) )-( (OMEGAN*XN) /  (  2 
OMEG=ABSF(OMEGAN)  - 

IF(0MEG-100o  )  15, 16, 16 

15  DMP  =  ABSF ( DAMPR  ) 

I F ( DM P-1 • ) 17, 16, 16 

16  WRITE  OUTPUT  TAPE  2,18 

18  FORMAT  ( 1 9H  SOLUTION  DIVERGING) 

DAMPR  =  .4  _  •  i- 

20  ALPHA=CAMPR*XDOO/XO 
BETA=SCRTF ( ALPHA**2.-XDDC0/X0 ) 

OMEGAN=-ALPHA-BETA _  _ 

weGPi0SRS«E^NrEG1N‘XD””'1  '°«6*S«^ar.-x6«n 

IFI0MEG-100. )21,22,?2  + - - - 

21  DMP=ABSF ( DAMPR ) 

IFIDMP-l. )23,22,22 
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22  WRITE  OUTPUT  TAPE  2t24  ..  _ 

24  FORMAT  (19H  SOLUTION  DIVERGING) 

GO  TO  12  _ 1 _ : _ 

23  WRITE  OUTPUT  TAPE  2,25,  OMEGAN,  DAMPR 

25  FORMAT  (  2F1  3.6  )  _ _ _ _ 

IF(J-20)26,12, 12 

26  J  =  J  +  T  _ _ _ 

GO  TO  20 

17  WRITE  OUTPUT  TAPE _2 ,  J 0.,  _0 M E G AN ,__D A  MP  R _ 

10  FORMAT  (2F13;6) 

I F ( J-2C ) 1 1,12,12 _ 

1 1  J  =  J+1  .... 

GO  TO  9  _ = _ 

12  J  =  1  . 

x  0 = x  n _ : _ 

X000=XCN 

X€DOO=XDDN  _ 

IFTK-4Q) 13,14, 14 

13  K=K+ 1  _ 

GO  TO  2 

14  CALL  EXIT  _ _ _ 

END  (  1  , 0,0, 0,0,0,  G,  0,0,  C  ,  C ,i>»Q»U»v) 
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